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FRAZIL  ICE  FORMATION 

R.  Ettema,  M.F.  Karim  and  J.F.  Kennedy 


INTRODUCTION 

Background 

Frazil  ice  forms  in  flowing  or  turbulent  water 
that  has  become  supercooled  by  heat  transfer  to 
overlying  air.  Turbulence,  by  agitating  the  water 
surface  or  the  water-air  interface,  inhibits  the 
growth  of  a  monolithic  ice  sheet.  The  principal 


phases  of  frazil  ice  growth  are  depicted  schemati¬ 
cally  in  Figure  1 .  Ice  crystals  such  as  snow,  frozen 
fog  or  ice  fragments  can,  on  entering  supercooled 
water,  initiate  the  formation  of  frazil  ice.  As  the 
crystals  of  frazil  ice  form,  they  are  entrained  and 
transported  downstream  in  flowing  water,  or  are 
dispersed  over  the  depth  of  agitated  water.  Once 
formed,  small  frazil  ice  crystals  grow,  agglutinate 


AGGLUTINATION,  SINTERING  AGGLOMERATION  OF  FRAZIL  AGGLOMERATION  OF  FLOCS 
NUCLEATION  OF  FRAZIL-ICE  DISKS  ICE  TO  FLOCS  TO  PANS  AND  SLUSH 


Figure  1.  Phases  of  frazil  ice  formation  in  a  stream  (after  Williams  1959). 
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to  each  other  and  may  subsequently  agglomerate  as 
floes.  The  buoyant  floes  of  frazil  ice  rise  to  the 
water  surface  where,  if  there  is  no  ice  cover,  they 
accumulate  to  form  ice  pans  on  their  exposure  to 
the  frigid  air.  The  pans  often  freeze  to  each  other 
and  form  small  ice  floes.  The  accumulation  of  these 
floes  along  the  edge  of  a  downstream  ice  cover  can 
lead  to  the  freeze-over  of  a  river.  When  the  trans¬ 
ported  floes  of  frazil  ice  become  trapped  beneath  an 
existing  ice  cover,  slush  ice  may  accumulate.  The  vol¬ 
ume  of  frazil  ice  produced  in  a  stream  can  potential¬ 
ly  exceed  the  volume  of  surface  ice  produced  in  it. 

Frazil  ice  is  of  engineering  significance  largely  be¬ 
cause  accumulations  of  it  can  cause  problems  for 
many  structures  in  or  adjoining  rivers,  lakes  and 
coastal  waters  in  cold  regions.  Two  prominent 
problems  associated  with  frazil  ice  are  the  adherence 
of  frazil  ice  to  structures  in  water  and  the  constric¬ 
tion  of  a  stream  or  river  flow  cross  section  by  accu¬ 
mulations  of  frazil  ice  beneath  an  ice  cover.  The 
present  study  directly  addresses  neither  of  these 
problems;  instead,  it  deals  with  some  of  the  funda¬ 
mental  aspects  of  frazil  ice  formation. 

Scope  of  study 

The  amount  of  frazil  ice  that  forms  in  a  body  of 
supercooled  water  principally  depends  on  the  degree 
of  supercooling  at  the  instant  of  seeding,  nucleation 
dynamics  and  the  rate  of  heat  transfer  from  the 
water  body.  In  nature,  these  factors  are  functions 
of,  among  other  things,  the  turbulence  characteris¬ 
tics  of  the  flow.  The  formation  of  frazil  ice  involves 
processes  of  crystal  growth,  secondary  nucleation  or 
the  creation  of  new  crystals  via  collision  breeding, 
and  the  combination  of  frazil  ice  crystals  to  form 
frazil  ice  particles  and  floes.  Three  critical  unre¬ 
solved  issues  related  to  the  formation  of  frazil  ice 
are: 

1 .  How  does  water  temperature  at  seeding  influ¬ 
ence  the  rate  of  frazil  ice  formation? 

2.  How  do  the  turbulence  characteristics  of  a 
flow  affect  the  rate  of  secondary  nucleation 
and  frazil  ice  formation? 

3.  How  does  the  mean  size  of  frazil  ice  particles 
vary  v/ith  water  temperature  and  the  turbu¬ 
lence  characteristics  of  the  body  of  water? 

The  study  reported  here  investigated  these  issues 
by  observation  of  them,  both  experimentally  in  a 
turbulence  jar  and  by  way  of  an  analytical  model  in 
which  the  rate  of  frazil  ice  formation  is  related  to 
turbulence  level  and  temperature  of  a  body  of  pure 
water.  The  results  from  the  experiments  guided  the 
formulation  of  the  analytical  model  that  predicts  the 
rate  of  frazil  ice  formation  for  a  given  temperature 


of  supercooling  and  turbulence  exchange  coefficient 
of  a  flow.  River  and  stream  waters  rarely  become 
supercooled  to  temperatures  less  than  about  -0.1  °C. 
However,  so  that  the  influences  of  turbulence  and 
temperature  of  a  flow  would  be  accentuated  suffi¬ 
ciently  to  be  apparent,  frazil  ice  formation  over  a 
wider  range  of  water  temperatures,  as  low  as  -0.9°C, 
was  studied. 

The  Literature  Review  section  of  this  report 
briefly  reviews  literature  pertaining  to  nucleation, 
growth  and  evolution  of  frazil  ice  particles  in  pure 
water.  Formulation  of  an  analytical  description  of 
frazil  ice  formation  is  outlined  in  the  Analytical 
Model  section.  Experiments  that  are  companion  to 
the  analytical  model  are  described  in  the  Experimen¬ 
tation  section.  Results  from  the  analytical  model 
and  from  experimentation  are  combined  and  dis¬ 
cussed  in  the  Results  section  and  the  principal  con¬ 
clusions  of  the  study  are  presented  in  the  Conclusions 
section.  Appendix  A  is  a  description  of  frazil  ice  ex¬ 
periments  that  were  performed  in  the  early  stage  of 
the  present  investigation,  but  were  subsequently 
abandoned.  A  computer  program  for  solving  the 
analytical  model  is  listed  in  Appendix  B  and  the  re¬ 
sults  obtained  from  this  analytical  model  are  pre¬ 
sented  in  Appendix  C. 

LITERATURE  REVIEW 
Introduction 

The  formation  of  frazil  ice  and  its  growth  and 
evolution  into  larger  ice  masses  in  rivers,  lakes, 
coastal  waters  and  man-made  channels  pose  many 
challenging  problems  for  the  design  and  operation 
of  hydraulic  structures.  Understanding  the  mechan¬ 
ics  of  frazil  ice  formation  has  been  the  subject  of  in¬ 
vestigation  by  many  researchers— some  of  the  earliest 
research  efforts  were  undertaken  by  Barnes  (1928) 
and  Devik  (1931).  Altberg  (1936)  presented  a  com¬ 
prehensive  review  of  early  work  (1915-1935)  on 
underwater  ice  formation  and  discussed  the  forma¬ 
tion  of  frazil  ice  and  anchor  ice  in  terms  of  heat 
convection  by  supercooled  water  (anchor  ice  is  sub¬ 
merged  ice  that  is  attached  to  the  bottom  of  a  stream 
or  submerged  objects  in  a  stream).  Tsang  (1982), 
Martin  (1981)  and  Osterkamp  (1978a)  have  recent¬ 
ly  published  erudite  reviews  of  current  knowledge 
concerning  frazil  ice.  The  present  review  is  limited 
to  surveying  briefly  the  literature  dealing  with  the 
salient  features  of  frazil  ice  formation,  namely,  the 
influences  of  water  temperature  and  turbulence  on 
the  rate  of  frazil  ice  formation  and  the  size  of  frazil 
ice  particles. 
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Incipient  formation  of  frazil  ice 

Frazil  ice  forms  in  supercooled  liquid  water  that 
is  turbulent.  The  introduction  of  seed  ice  crystals  is 
required  to  begin,  by  the  process  of  secondary  nucle- 
ation,  the  formation  of  frazil  ice.  Secondary  nuclea- 
tion  involves  the  proliferation  of  ice  crystals  through 
the  mechanism  of  collision  and  disintegration  of  an 
initial  set  of  ice  particles.  Ice  particles  may  break 
during  collisions  with  each  other  or  with  other  solid 
surfaces,  or  because  of  their  impact  with  fluid  eddies. 

Depending  on  the  rate  of  heat  transfer  to  overly¬ 
ing  air,  a  body  of  water  may  be  further  supercooled 
after  frazil  ice  has  begun  to  form,  or  it  may  begin  to 
warm  because  of  the  release  of  latent  heat  of  fusion 
that  accompanies  ice  formation.  In  nature,  the  super¬ 
cooling  of  water  is  usually  accompanied  by  ice 
growth,  either  as  frazil  ice  or  an  ice  cover  forming 
from  the  surface  borders  of  the  flow.  Consequently, 
the  supercooling  of  water  in  rivers  and  streams  is 
usually  limited  to  about  0.1  °C  or  less  below  0°C. 
Williams  (1959),  in  his  review  of  the  current  theories 
of  frazil  ice  formation,  summarized  various  research¬ 
ers’  theories  of  frazil  ice  formation.  These  studies 
dealt,  to  varying  degrees,  with  the  supercooling  of 
water  and  the  nucleation  of  supercooled  water,  and 
included  studies  concerning  the  factors  affecting  the 
rate  of  cooling  in  open  rivers.  In  addition,  early  at¬ 
tempts  to  predict  the  onset  and  rate  of  frazil  ice  pro¬ 
duction,  based  on  knowledge  of  water  temperature, 
are  outlined  by  Williams.  The  following  citations 
paraphrase  a  portion  of  Williams’  paper. 

Schaefer  (1950)  reported  that  the  supercooling 
of  water  in  streams  frequently  persisted  for  long 
periods  (often  days)  but  that  the  water  is  rarely  super¬ 
cooled  below  -0.1°C.  According  to  Dorsey  (1948) 
the  formation  of  frazil  ice  is  started  by  ice  molecules 
or  by  heterogeneous  singularities— that  is,  foreign  par¬ 
ticles  in  water  that  serve  as  ice  nuclei-that  fall  into  a 
supercooled  body  of  water.  Kumai  and  Itagaki  (1954) 
concluded  from  their  photographic  study  of  frazil  ice 
formation  that  with  a  temperature  of  supercooling 
down  to  -0.3°C  discoid  and  spicule  frazil  ice  particles 
are  produced,  and  from  -0.6°  to  -0.9°C  only  the 
spicule  particles  are  formed.  Granbois  (1953)  report¬ 
ed,  from  his  measurements  in  rivers,  that  frazil  ice 
forms  only  when  the  rate  of  temperature  change  of 
the  water  is  greater  than  0.01°C/hr  for  temperatures 
of  supercooling  of  -0.1  °C  or  higher— an  ice  sheet 
forms  over  the  water  surface  for  lower  rates  of  tem¬ 
perature  change. 

Wigle  (1970)  and  Arden  and  Wigle  (1972)  describe 
a  study  in  the  Niagara  River  where  frazil  ice  produc¬ 
tion  began  at  a  water  temperature  of  -0.05°C  during 
cold,  clear  nights,  and  developed  into  “a  driving  snow¬ 
storm  such  as  can  be  seen  through  the  headlight  beams 


of  a  moving  automobile  at  night.”  They  estimated 
that  the  accumulated  frazil  ice  reduced  the  cross- 
sectional  areas  of  the  river  by  as  much  as  25%.  Ex¬ 
periments  by  Garabedian  and  Strickland-Constable 
(1974)  shed  light  on  the  rapid  multiplication  of  frazil 
ice  nuclei  through  the  mechanism  of  “collision  breed¬ 
ing”;  at  a  supercooling  of  -0.1 0°C,  44  crystals  formed 
in  5  seconds  from  one  original  seed  crystal. 

Osterkamp  (1978a)  observed  that  the  nucleation  of 
frazil  ice  in  rivers  is  heterogeneous  for  small  values  of 
temperature  difference  below  the  freezing  point  of 
water  ATn  (or  supercooling).  In  addition,  he  con¬ 
cluded  that  flow  turbulence  per  se  does  not  notice¬ 
ably  affect  ATn  and  that  water  of  an  ice-free  river 
does  not  supercool  to  less  than  -0.05°C  at  depths  of 
10-mm  and  more  below  the  water  surface.  In  an¬ 
other  paper,  Osterkamp  (1978b)  reviewed  various 
mechanisms  for  frazil  ice  nucleation  and  concluded 
that  mass  exchange  and  secondary  nucleation  are  the 
most  likely  mechanisms.  The  former  mechanism  in¬ 
volves  ice  particles  that  originate  at  some  distance 
from  the  water  body  (snow,  frost,  ice  particles  from 
banks,  trees,  etc.)  and  ice  particles  that  originate  in 
the  river  water;  for  example,  water  droplets  or  vapor 
can,  on  being  propelled  into  the  cold  air  above  the 
water  surface  by  splashing,  wind  spray  and  evapora¬ 
tion,  freeze  and  fall  back  to  the  water  surface  after 
being  spontaneously  nucleated  in  the  cold  air. 

Mueller  (1978)  conducted  several  experiments 
with  temperatures  of  supercooling  between  -0.05° 
and  -0.30°  C  in  a  turbulence  jar  (17.2  by  12  by  20 
cm),  which  was  an  earlier  version  of  the  jar  in  which 
the  present  experiments  were  performed.  Water  in 
the  jar  was  cooled  below  the  freezing  point  by  circu¬ 
lation  of  a  coolant  through  tubes  in  the  walls  and 
base  of  the  jar  and  turbulence  was  generated  by  a 
moving  grid.  In  Mueller’s  experiments  no  frazil  ice 
formed  in  the  turbulence  jar  unless  the  supercooled 
water  was  seeded  with  ice  nuclei.  He  found  that  this 
happened  regardless  of  the  degree  of  supercooling, 
the  intensity  of  turbulence  or  the  presence  of  foreign 
materials  in  the  water. 

In  summary,  the  weight  of  both  laboratory  and 
field  evidence  strongly  indicates  that  the  heterogene¬ 
ous  nucleation  of  supercooled  water  by  the  ice  frag¬ 
ments  or  crystals  initiates  the  formation  of  frazil  ice. 
Nevertheless,  Hobbs  (1974)  alludes  to  the  nucleation 
of  supercooled  water  by  the  immersion  of  a  cold 
metal  rod  in  the  water. 

Particle  size  and  evolution  of  frazil  ice 

Notwithstanding  the  many,  although  somewhat 
diffuse,  descriptions  of  size  and  shape  of  frazil  ice 
particles,  the  correlation  of  frazil  ice  geometry  with 
water  temperature  has  not  been  well  delineated.  In 


addition,  a  rather  varied  terminology  is  used  in  pub¬ 
lished  literature  to  describe  frazil  ice  particles. 

Field  observations  of  frazil  ice  growth  in  rivers, 
as  reported  by  Arakawa  (1954),  Michel  (1967), 

Arden  and  Wigle  (1972),  Osterkamp  (1978b)  and 
others,  indicate  that  frazil  ice  formation  begins  at 
a  temperature  of  supercooling  of  the  order  of  -0.05° 
to  -0.1°C,  with  the  frazil  ice  having  the  form  of  cir¬ 
cular  disks,  1  to  5  mm  in  diameter.  Schaefer  (1950) 
recorded  the  dimensions  of  frazil  ice  particles  in  the 
Mohawk  River.  He  found  that  disk  diameter  ranged 
from  1  to  5  mm  and  disk  thickness  varied  from 
2.5  x  10'2  to  10"'  mm.  Hobbs  (1974)  concluded 
that  the  circular  disk  is  the  preferred  growth  form 
for  frazil  ice  forming  in  pure  water  that  has  been 
supercooled  to  within  a  temperature  range  of  -0.1° 
to  -0.3°C.  The  concensus  of  those  who  have  ob¬ 
served  frazil  ice  growth  appears  to  confirm  Hobbs’ 
conclusion. 

Laboratory  experiments  by  Carstens  ( 1 966), 
Hobbs  (1974),  Mueller  (1978)  and  other  research¬ 
ers  report  particle  sizes  for  frazil  ice  ranging  up  to 
10  mm  in  diameter.  The  larger  sizes  reported  for 
frazil  ice  particles  grown  in  laboratories  are  probably 
attributable  to  the  lower  water  temperatures  that 
can  be  attained,  and  used,  in  laboratory  experiments. 
This  observation,  however,  remains  to  be  verified. 

The  evolution  of  frazil  ice  into  larger  agglomer¬ 
ates,  according  to  Michel  (1967, 1972)  and  Oster¬ 
kamp  (1978a),  starts  when  crystals  of  1  to  5  mm  in 
diameter  and  10_2to  KT1  mm  in  thickness  form 
“floes”  or  collections  of  disks  5  to  100  mm  in  diam¬ 
eter;  subsequendy  “pans”  with  diameters  of  the 
order  of  1  m  and  thickness  of  0.1  to  0.5  m  form; 
“floes”  with  diameters  of  1  to  30  m  and  thickness 
of  0.5  to  5  m  may  then  form  as  patchworks  of  pans. 
Carstens’  (1970)  analysis  of  several  Norwegian  rivers 
showed  that  for  surface  velocity  vs  greater  than  0.6 
m/s,  but  less  than  1.2  m/s,  frazil  ice  forms  and  accu¬ 
mulates  on  the  water  surface-for  vs  greater  than 
1 .2  m/s  frazil  ice  forms  over  the  whole  river  depth. 
These  values  are  roughly  in  agreement  with  experi¬ 
ments  conducted  by  Hanley  and  Michel  (1977),  who 
also  found  0.24  m/s  to  be  the  threshold  velocity  for 
frazil  ice  formation. 

Several  anecdotal  accounts  of  frazil  ice  accumula¬ 
tion  in  rivers  and  lakes  have  been  published.  A  1 6- 
km-long  under-ice  accumulation  of  frazil  ice  with  an 
approximate  volume  of  5  x  107  m3  in  the  La  Grande 
River  is  described  by  Michel  (1978).  He  observed 
that  the  accumulation  began  with  the  collection  of 
frazil  ice  beneath  the  solid  ice  cover,  which  was 
some  distance  downstream  from  rapids.  The  open 
water  at  the  rapids  enabled  the  water  to  become 
supercooled.  Beginning  at  a  downstream  location, 


the  accumulation  then  progressed  upstream.  Michel 
(1978)  and  Tesaker  (1975)  developed  a  Froude  num¬ 
ber  criterion  to  predict  the  maximum  thickness  of 
frazil  deposition  beneath  the  ice  cover. 

Influences  of  turbulence  and  water  temperature  on 
the  rate  of  frazil  ice  formation 

Although  a  number  of  analytical  models  of  frazil 
ice  growth  have  been  formulated  and  several  experi¬ 
mental  studies  of  frazil  ice  growth  conducted,  no 
study  has  adequately  explained  the  influence  of  tur¬ 
bulence  and  water  temperature  on  the  rate  of  frazil 
ice  growth. 

One  of  the  earliest  models  of  frazil  ice  growth 
was  proposed  by  Carstens  (1966),  who  suggested  the 
following  relationship  for  the  rate  of  frazil  ice  growth, 
dMJdt  (Afj  =  mass  of  ice,  t  =  time): 

dM-  j  dT 


where  Lt  =  latent  heat  of  fusion 
5h  =  specific  heat  of  water 
T  -  water  temperature 
Mw  =  mass  of  water 

Qs  =  rate  of  heat  loss  from  water  surface. 

Another  expression  for  dMJdt,  from  the  point  of 
view  of  nucleation  and  growth  of  ice  crystals,  was 
proffered  by  Carstens  for  the  rate  of  frazil  ice  pro¬ 
duction: 

dM- 

L^=aA^Tm-T)  (2) 

where  A  t  =  surface  area  of  crystals  at  time  t ,  7^  = 
freezing  temperature  of  water,  and  a  =  heat  transfer 
coefficient,  given  by 


k 

0t=JN'1 


where  k  =  thermal  conductivity  of  water,  C  =  charac¬ 
teristic  dimension  of  ice  crystals,  and  Nu  =  Nusselt 
number,  which  governs  the  heat  transfer  from  ice 
surfaces.  /Vu  depends  on  the  geometry  of  ice  crystals 
and  on  the  structure  of  flow  around  the  ice  crystals. 
For  a  given  geometry  of  ice  crystal,  Carstens  proposed 
that 


NU=C{R? 


where  C,  =  constant 


Re  =  Reynolds  number  =  vi/v 


;$ 

H 


v  =  a  characteristic  velocity  of  the  flow 
v  =  kinematic  viscosity 
m  -  constant  (0.5  ~~  0.8). 

Carstens  further  proposed  that 


At  ~  Al/3  M2'3 

and 
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C  ~  ,\~113  Allp  (6) 

where  (V  =  number  of  ice  crystals.  By  utilizing  eq 
3  through  6,  we  can  rewrite  eq  2  as 

y-  A/V3  A,2/3/R"'(7m-7')  ■  (7) 

Carstens  concluded  from  observations  of  ice  crystal 
formation  in  his  experimental  flume  that  the  number 
of  ice  crystals  N  increases  with  increasing  rate  of 
heat  loss  Qs.  This  observation  was  in  agreement 
with  his  analytical  formulation,  expressed  in  eq  1 
through  7.  An  important  observation  indicated  by 
his  experimental  and  analytical  data  is  that,  for  de¬ 
creasing  values  of  the  Reynolds  number,  the  re¬ 
quired  condition  for  prolonged  periods  of  supercool¬ 
ing  can  be  determined  from  the  relationship 

R”‘  IV213  Mj13  =  constant.  (8) 

In  addition  to  formulating  the  above  analytical 
model  of  frazil  ice  formation,  Carstens  also  conduct¬ 
ed  experiments  in  a  race-track  shaped  flume  (1.2  by 
0.3  m  by  0.2  m  wide)  at  the  Technical  University  of 
Norway.  The  bottom  and  sides  of  the  flume  were 
insulated,  thereby  allowing  heat  loss  only  through 
the  free  surface  to  the  air  stream  above.  The  air 
temperature  was  maintained  at  -10°C.  With  a  cool¬ 
ing  rate  of  0.01°C/minute,  he  attained  a  maximum 
supercooling  of  -0.04°C.  The  rate  of  cooling  was 
varied  by  changing  the  temperature  and  the  distance 
between  the  flume  and  a  fan  that  controlled  air 
speed.  His  experimental  results  show  that  an  in¬ 
crease  in  the  rate  of  heat  loss  increased  the  tempera¬ 
ture  of  supercooling  attained,  the  residual  supercool¬ 
ing  temperature  and  the  rate  of  temperature  rise 
from  maximum  supercooling. 

The  results  of  Mueller’s  (1978)  experiments  show 
that  the  growth  rates  of  frazil  ice  particles  increase 
with  greater  supercooling  and  higher  grid  velocity 
(that  is,  with  increased  turbulence).  From  his  ob¬ 
servations  that  the  heat  transfer  per  particle  could 
be  normalized  with  supercooling  and  that  the  size  of 
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particles  was  constant  in  all  experiments,  Mueller 
concluded  that  the  rate  and  the  total  amount  of 
frazil  ice  production  can  be  predicted  if  estimates 
of  tiie  number  of  particles  and  heat  transfer  per 
particle  can  be  made.  He  suggested  that  the  normal¬ 
ized  heat  transfer  can  be  related  to  the  turbulent 
characteristic  of  flow,  but  a  nucleation  theory  was 
not  avauaole  to  estimate  the  number  and  growth 
rate  of  frazil  particles.  His  experiments  also  indicat¬ 
ed  that  frazil  production  increases  with  increasing 
grid  frequency  or  level  of  turbulence  intensity. 

Andres  (1982)  carried  out  an  experimental  pro¬ 
gram  to  study  the  effects  of  turbulence  and  the  rate 
of  heat  loss  on  the  supercooling  of  water  and  the 
subsequent  frazil  ice  formation.  He  performed  his 
experiments  using  an  octagonal  tank  at  the  base  of 
which  was  mounted  a  propeller  to  generate  a  cyclo- 
strophic  flow  in  which  the  frazil  ice  formed.  From 
his  data,  Andres  observed  a  trend  that  indicated 
that,  for  increasing  values  of  turbulence  intensity, 
the  period  of  supercooling  decreased  but  the  rate  of 
frazil  ice  formation  increased.  The  roles  played  by 
turbulence  and  water  temperature  in  the  rate  of 
frazil  ice  formation  are  difficult  to  define  from 
Andres’  data. 

An  elaborate  mathematical  model  of  frazil  ice 
formation  and  transport  in  a  river  with  a  steady  flow 
was  constructed  by  Matousek  (1981 ).  His  model  is 
formulated  using  a  heat  balance  equation  for  a  steady 
uniform  flow.  In  his  state-of-the-art  review  of  re¬ 
search  concerning  frazil  ice,  Tsang  (1982)  derived 
the  following  relationship  for  the  rate  of  ice  produc¬ 
tion: 


=  £h  Q,-Qi 

dr  /.j  dr  1  pip 


where 


fj  =  concentration  of  frazil  by  weight 
SH  =  specific  heat  of  water 
T  =  water  temperature 
Z-j  =  latent  heat  of  fusion  of  ice 
Q%  =  heat  flux  from  water  to  air 
Oj  =  heat  flux  from  channel  bed  to  water 
D  =  flow  depth. 


By  noting  that  dcjdt  =0just  prior  to  nucleation, 
we  may  rewrite  eq  9  as 


dCj 

dr 


h 


(dt  n 


(10) 


where  ( dT/dr)n  =  dT/dt  at  nucleation.  Integration 
of  eq  10  yields 


SH  7\T 


ANALYTICAL  MODEL 


where  Tn  =  nucleation  temperature  and  rn  =  the 
time  of  first  nucleation.  Tsang  comments  that,  as 
seen  from  eq  1 1,  frazil  concentration  (Cj)  is  greatly 
influenced  by  (3 T/dt)n  which,  in  turn,  is  deter¬ 
mined  by  heat  fluxes  to  and  from  the  water  body. 

Analytical  models  such  as  those  constructed  by 
Carstens  (1966),  Matousek  (1981)  and  Tsang  (1982) 
are  useful  in  that  they  may  indicate  the  nature  of 
frazil  ice  formation  and  its  relationship  with  water 
temperature  and  turbulence.  However,  a  drawback 
of  these  models  is  that  they  are  rather  unwieldy  and 
difficult  to  use. 

Conclusions 

Despite  the  burgeoning  amount  of  work  that  has 
been  invested  in  the  study  of  frazil  ice,  the  three 
issues  outlined  in  the  Scope  of  Study  section  have 
remained  largely  unaddressed.  These  are  the  influ¬ 
ences  of  both  water  temperature  at  nucleation  and 
turbulence  characteristics  of  a  flow  on  the  rate  of 
frazil  ice  formation  and  the  variation  of  frazil  ice 
particle  size  with  the  temperature  and  turbulence  of 
a  flow. 

Evidence  from  reported  observations  suggests 
that  natural  rivers  and  streams  rarely  become  super¬ 
cooled  to  temperatures  less  than  about  -0.1°C. 

Frazil  ice  particles  formed  at  this  temperature  tend 
to  be  disks  with  diameters  typically  attaining  a  max¬ 
imum  value  of  approximately  2  to  3  mm,  and  thick¬ 
nesses  varying  from  1(T2  to  KT1  mm.  Larger  diam¬ 
eters  of  frazil  ice  particles  are  reported  from  labora¬ 
tory  studies,  frequently  involving  lower  water  tem¬ 
peratures  than  are  recorded  for  rivers  and  streams, 
thus  suggesting  that  particle  size  is  a  function  of 
water  temperature.  The  influence  of  turbulence  on 
particle  size  appears  not  to  have  been  studied.  The 
amount,  or  concentration,  of  frazil  ice  that  forms  in 
a  flow  is  indicated  by  several  researchers  to  increase 
with  turbulence  intensity  of  a  flow,  but  it  is  unclear 
if  this  occurs  by  turbulence  causing  the  water  to  be¬ 
come  supercooled  to  lower  temperatures  as  well  as 
by  playing  a  direct  role  in  the  physics  of  frazil  ice 
formation.  Field  data  (for  example  those  reported 
by  Tsang  [1982] )  indicate  that  an  upper  limit  to 
the  concentration  of  frazil  ice  that  can  be  formed  in 
supercooled  water  of  streams,  rivers  and  lakes  is  of 
the  order  of  0.5%  by  weight. 


Introduction 

By  combining  elements  of  heat  transfer  theory 
associated  with  the  formation  of  ice  in  liquid  water 
with  elements  of  turbulence  theory,  as  it  applies  in 
a  turbulence  jar,  it  is  possible  to  identify  the  influ¬ 
ences  of  water  temperature  and  turbulence  on  the 
rate  of  frazil  ice  formation.  Experimentation  was 
used  to  gain  knowledge  of  the  influences  of  turbu¬ 
lence  and  water  temperature  on  the  temporal  varia¬ 
tion  of  frazil  ice  particle  size.  Turbulence  was  gener¬ 
ated  in  the  turbulence  jar  by  a  vertically  oscillating 
grid  that  had  to  be  heated  to  stop  frazil  ice  from  ad¬ 
hering  to  it. 

The  analytical  model  is  developed  from  a  Lagran- 
gian  standpoint-that  is,  the  formation  of  frazii  ice 
is  considered  for  a  supercooled  parcel  or  volume  of 
water  moving  with  the  flow  in  a  stream,  river  or  lake. 
The  analysis  deals  with  the  rates  and  the  amounts  of 
frazil  ice  that  can  form  in  a  parcel  of  water  that  has 
been  supercooled  to  a  temperature  T0  prior  to  the 
formation  of  frazil  ice. 

Elements  of  heat  transfer 

When  ice  is  formed  in  supercooled  water,  latent 
heat  of  fusion  is  released,  thereby  raising  the  tem¬ 
perature  of  the  surrounding  water.  The  rate  of 
frazil  ice  formation  can  be  related  to  the  rate  of 
latent  heat  release  in  order  to  provide  a  framework 
for  the  study  of  the  influences  of  water  temperature 
and  turbulence  on  the  rate  of  frazil  ice  formation. 
Formulation  of  the  analytical  framework  proceeds 
as  follows.  The  rate  at  which  heat  is  taken  up  by 
the  water  is  given  by 

Qw=PSHVj-t(T-TJ  (12) 

where  p  =  water  density 

SH  =  specific  heat  of  water 
V  -  volume  of  water 
T  =  water  temperature 
Tm  =  equilibrium  temperature  of  liquid 

water  and  ice  (for  pure  water,  Tm  =0°C). 

The  rate  at  which  heat  is  transferred  to  the  water 
from  its  surroundings  and  from  the  heated  agitator 
grid  is 

Qh=pSHVT  (13) 


6 


where  T  =  the  rate  of  temperature  rise  attributable 
to  these  heat  sources.  It  was  experimentally  found 
that  T  could  be  taken  to  be  constant  for  a  given 
combination  of  ambient  air  temperature  and  grid 
temperature.  The  temperature  rise  terms  are  illus¬ 
trated  in  Figure  2. 


where*  =  (7"-  rm), and 


Q  r  pSHV 


_“ZL  .  r=  a 
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WATER  TEMP 
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For  pure  water  and  the  Celsius  temperature  scale 
(Tm  =  0),  *  =  T,  and  thus  eq  1 6  becomes 


=  T  -(Bt  + A)  T- CT 2 . 


Figure  2.  Temporal  variation  of  temperature  in  the 
turbulence  jar  (Ta  =  water  temperature  at  seeding). 

The  surface  area  of  frazil  ice  crystals,  At,  at  any 
time  t  after  seeding  can  be  expressed  as 

pSHV 

A'  =A°  +~p^e  (r"r°  ‘ 773  (14) 

where  Aa  =  surface  area  of  ice  at  t  =  0  (the  time  of 
seeding) 

Pi  =  density  of  ice 
Zj  =  latent  heat  of  fusion  of  ice 
8  =  characteristic  thickness  scale  of 
the  ice  particles 

Ta  =  supercooling  at  the  time  of  seeding. 

As  frazil  ice  particles  are  flat  platelets,  the  surface 
area  of  ice  need  not  be  reckoned  in  terms  of  particle 
shape  and  number  of  particles.  Rather,  the  particles 
can  be  considered  to  be  fragments  of  a  single,  thin 
ice  sheet,  with  heat  being  transfered  from  both  faces 
of  the  sheet.  Heat  transfer  from  the  edges  of  the 
sheet  can  be  neglected.  The  surface  area  in  eq  14 
represents  the  area  of  both  faces  of  the  ice  sheet. 

The  rate  at  which  heat  of  fusion  is  given  up  by  the 
ice  particles  is 

pS„V 

<?f.«M0+^<r-r0-n>i(rm-n 

(15) 

where  a  =  coefficient  of  heat  transfer  between  the 
frazil  particles  and  water. 

Equating  to  the  sum  of  @h  and  Qf  yields 


-=  T  -  {Bt  +  A) x  -  Cx2 
at 


If  the  surface  area  of  the  edges  of  a  frazil  ice 
platelet  (the  platelets  typically  are  about  0.05  to 
0.10  mm  thick)  is  neglected,  the  characteristic  width 
C  can  be  taken  to  be  b/ 2  and  the  frazil  ice  concentra¬ 
tion  at  time  t,  ct,  is  expressed  as 

A  fb  Anb  pS\i  • 

c'  =  tv--fv+^T-T°-Tt)  (I8) 

where  b  =  representative  thickness  of  frazil  ice 
platelets,  which  was  assumed  to  remain  constant  at 
a  value  of  0.1  mm  throughout  each  experiment. 

Note  that  T  can  be  removed  from  eq  17  and  18 
to  yield  a  differential  equation  for  frazil  ice  concen¬ 
tration  ct .  It  is  clear,  both  from  physical  reasoning 
and  eq  18,  that  c,  varies  linearly  with  (T-  Tn  -  ft), 
or  if  T  and  T  are  both  set  to  zero,  c,  varies  linearly 
with  temperature  of  supercooling  T0. 

Several  refinements  to  the  above  analysis  can  be 
considered.  One  refinement  is  to  allow  the  heat 
transfer  coefficient  a  to  decrease  with  increasing 
concentration  of  frazil  ice.  As  the  concentration  of 
frazil  ice  increases,  the  effective  dynamic  viscosity 
of  the  water-ice  mixture  increases  because  the  super¬ 
cooled  water  that  was  initially  clear  ultimately  be¬ 
comes  a  slush.  The  following  heuristic  relationship 
was  formulated  to  describe  the  variation  of  a  with 
the  concentration  of  frazil  ice  ct : 

a  =  a0n,  (1  +n2  e  (”3C^)  (]9) 

where  aQ  =  initial  value  assumed  for  the  heat  trans¬ 
fer  coefficient  when  concentration  ct  is  zero,  and 
nt ,  n2  and  n3  are  experimentally  determined  con¬ 
stants. 

Equation  1 7  was  solved  for  several  values  of  aQ 
in  order  to  give  the  temporal  values  of  mean  water 
temperature  and  concentration  of  frazil  ice  for 
supercooled  water.  To  facilitate  the  solution,  A0 
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was  assumed  to  be  equivalent  to  a  single  frazil  ice 
crystal  with  two  0.25-cm2  surfaces  and  a  thickness 
of  0.1  mm.  This  is  typical  of  the  size  of  frazil  ice 
particles  that  were  observed  during  the  experiments, 
and  gives  a  volume  equivalent  to  the  piece  of  seeding 
ice.  Results  were  produced  from  the  analytical 
model  by  way  of  a  computer-based  Runge-Kutta 
solution  routine,  with  the  Gill  modification.  A 
listing  of  the  computer  program  used  to  facilitate 
the  solution  is  given  in  Appendix  B. 

The  foregoing  analysis,  eq  12  through  19,  is  a 
global  description  of  frazil  ice  formation  in  a  parcel 
of  supercooled  water.  It  does  not  account  for  the 
influences  of  the  temporal  variation  of  size  and  size 
distribution  of  frazil  ice  on  the  rate  and  amount  of 
it  formed  in  a  parcel  of  supercooled  water.  Addi¬ 
tional  limitations  of  the  analysis  include  the  lack  of 
detailed  relationships  among  the  growth  of  frazil 
ice  crystals,  secondary  nucleation  and  heat  transfer 
from  suspended  crystals  in  a  turbulent  fluid.  None¬ 
theless,  the  analysis  elucidates  the  roles  played  by 
water  temperature  and  turbulence  level  on  the  rate 
and  amount  of  frazil  ice  formed  in  a  body  of  super¬ 
cooled  water.  It  also  provides  a  suitable  framework 
into  which  further  analytical  refinements  can  be 
added. 

Elements  of  turbulence 

The  relationships  between  the  turbulence  struc¬ 
ture  of  a  flow  and  the  formation  of  frazil  ice  or, 
more  pertinently,  between  turbulence  and  heat 
transfer  from  a  frazil  ice  particle,  are  involute  and 
not  well  defined.  Essentially,  the  action  of  turbu¬ 
lence  is  to  increase  the  rate  of  heat  transfer  from 
frazil  ice  particles  and  thereby  to  increase  the  rate 
of  frazil  ice  formation.  This  happens  provided  that 
the  particles  are  considerably  larger  than  the 
Kolmogorov  length  scale  or  microscale  of  the  turbu¬ 
lence  field  in  which  they  form,  and  that  the  frazil 
ice  particles  are  well  spaced,  so  that  they  can  be 
treated  as  single  particles  in  a  shear  flow. 

The  heat  transfer  from  small  ice  particles  to  their 
surrounding  fluid  depends  on  the  ratio  of  the  micro¬ 
scale  of  turbulence  to  the  size  of  the  particle.  Ice 
particles  are  convected  with  the  fluid  without  sens¬ 
ing  turbulence  if  the  microscale  is  large  compared  to 
the  particle.  Consequently,  the  relative  motion  and 
the  heat  transfer  between  the  particle  and  fluid  are 
small.  As  the  ratio  of  turbulence  microscale  to 
particle  size  decreases,  the  rates  of  relative  motion 
and  heat  transfer  increase.  In  addition,  as  long  as 
the  level  of  turbulence  is  such  that  the  microscale 
for  time  is  short  compared  to  the  time  scale  of  heat 
conduction,  the  heat  transfer  from  a  frazil  ice  parti¬ 
cle  to  the  fluid  will  be  an  order  of  magnitude  higher 


than  if  the  fluid  were  at  rest.  The  relative  magnitudes 
of  convective  heat  transfer  and  heat  transfer  by  con¬ 
duction  are  usually  expressed  in  terms  of  the  Nusselt 
number  N„ , 


N  =2* 
u  k 


where  k  =  coefficient  of  thermal  conductivity  and 
L  =  a  characteristic  length.  Nusselt  numbers  are  fre¬ 
quently  related  to  Reynolds  and  Prandtl  numbers, 
Re  and  Pr ,  respectively;  for  example 


Hu=KRne  P? 


where  n  and  m  =  exponents 
K  =  constant 

SH  =  specific  heat  of  the  fluid 
u  =  kinematic  viscosity  of  the  fluid 
p  =  density  of  the  fluid. 

Equation  21  is  primarily  used  to  evaluate  the  heat 
transfer  coefficient  a  for  a  particular  flow  case. 
Carstensf  1966)  and  Mueller  (1978)  adopted  lines 
of  analysis  involving  eq  21  in  order  to  estimate  the 
rate  of  heat  transfer,  and  thus  the  rate  of  frazil  ice 
formation,  for  various  levels  of  turbulence  in  a  tur¬ 
bulence  jar.  Because  of  the  lack  of  a  well-developed 
theory  to  describe  the  motion  of,  and  heat  transfer 
from,  discoidal  particles  in  a  shear  flow,  the  lines  of 
analysis  involving  relationships  of  the  form  of  eq 
21  as  followed  both  Carstens  and  by  Mueller  are 
fraught  with  many  difficulties. 

The  approach  followed  in  the  analysis  conducted 
for  this  study  was  to  use  the  principles  of  particle 
suspension  in  a  turbulent  flow  to  obtain  estimates 
of  turbulence  level,  then  to  obtain  values  of  heat 
transfer  coefficient  a  by  calibrating  the  analytical 
model  of  frazil  ice  formation  with  the  experimental¬ 
ly  determined  values  of  turbulence  exchange  coeffi¬ 
cient  e. 

Drawing  on  an  analogy  with  Prandtl’s  mixing- 
length  theory  (1925)  for  turbulence,  Rouse  (1937) 
described  the  suspension  of  particles  in  a  flow  in  the 
following  manner:  “If  at  any  flow  depth  there  are 
n  particles  per  unit  volume  of  fluid,  the  concentra¬ 
tion  gradient  being  written  [as]  dn/dy,  then  due  to 
transverse  fluctuations,  fluid  masses  bearing  n  parti¬ 
cles  per  unit  volume  will  be  carried  a  mean  distance 
[Cp  Prandtl  mixing  length]  across  the  flow  to  regions 
where  the  concentration  differs  by  the  amount 
[£p  dn/dy] .  Thus,  the  temporal  rate  of  passage  of 
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sediment  per  unit  area  [/Vs]  will  be  the  product  of 
the  rate  per  unit  area  of  transverse  flow  [|v'vl|  and 
the  difference  between  the  sediment  concentration 
in  the  traveling  fluid  mass  and  that  of  the  region  in¬ 
to  which  it  conies  (- C  Jn/Jy]."  Rouse  wrote 


The  variation  of  suspended  particle  concentration 
with  depth  can  be  related  to  the  turbulence  charac¬ 
teristics  of  the  How  by  using  the  following  relation¬ 
ships  proposed  by  Rouse  (1938):  “For  sediment 
suspension  it  is  customary  to  assume  that  a  state  of 
equilibrium  [exists)  between  the  mean  vertical  rate 
of  sediment  transport  due  to  turbulence  and  the  nor¬ 
mal  rate  of  settling  of  sediment  due  to  its  own  weight. 
Thus,  if  »■  denotes  the  settling  velocity  [of  the  sus¬ 
pended  particles  |  and  [t-  |  is  the  sediment  concen¬ 
tration  per  unit  volume  of  the  suspension 


Integration  leads  at  once  to  an  expression  for  the 
relative  concentration  [of  suspended  particles]  at 
any  point  above  some  arbitrary  reference  level  a, 


111  (  —  )  =  -  '  J  /  T7T 
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the  magnitude  of  [cM  |  depending  ultimately  upon 
conditions  in  the  boundary  region.”  In  the  foregoing 
quotations,  the  bracketed  portions  and  the  equation 
numbers  are  inserted  by  the  authors.  It  is  pertinent 
to  note  in  Rouse’s  analysis  of  particle  suspension  by 
turbulence  that  the  turbulence  momentum  exchange 
coefficient  e  is  equated  to  0es,  where  es  is  the  tur¬ 
bulence  mass  transfer  coefficient  and  (3  is  a  propor¬ 
tionality  coefficient,  which  for  fine  sediments  (fine 
sands  and  silts)  is  assumed  to  be  of  the  order  of  unity. 
For  coarser  sediments  0  may  decrease  with  increasing 
particle  size  and  no  longer  be  a  constant  value,  al¬ 
though  this  is  still  a  matter  of  some  controversy. 

The  turbulence  jar,  being  a  tank  in  which  a  grid 
oscillates  with  simple  harmonic  motion,  produces  a 
flow  with  essentially  the  same  degree  of  mixing  (v'2p) 
that  is  constant  over  the  major  portion  of  the  fluid 
volume.  The  turbulence  jar  is,  therefore,  well  suited 
for  the  study  of  heat  transfer  and  frazil  ice  formation 
in  a  turbulent  flow.  Integration  of  eq  24  for  constant 
v'8p  yields 


in  which  the  turbulence  exchange  coefficient  e  is 
equated  to  0  v'8p.  Rearranging  eq  25  leads  to 

(2« 
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from  whence,  noting  that  the  left-hand  side  of  eq 
26  is  equal  to  the  slope  of  the  concentration  versus 
depth  curve  S, 

e  =  (3v'fip  =  g.  (27) 

By  measuring  the  distribution  of  suspended  sediment 
concentration  over  the  depth  of  the  turbulence  jar, 
it  is  possible  to  relate  the  characteristic  magnitude  of 
the  turbulence  exchange  coefficient  e  to  the  frequen¬ 
cy  of  grid  oscillation  F  A  relation  between  e  and 
F  is  developed  from  experimental  results  in  the  Re¬ 
sults  section. 


EXPERIMENTATION 
Experimental  apparatus 

Several  designs  of  experimental  apparatus  to  pro¬ 
duce  a  turbulent  flow  field  were  investigated  during 
the  course  of  the  study.  The  experience  gained  from 
using  various  flow  devices  to  produce  a  turbulence 
field  f  or  study  of  the  effect  of  turbulence  and  water 
temperature  on  frazil  ice  formation  led  to  the  formu¬ 
lation  of  the  following  guidelines  for  experiments  of 
this  kind: 

1.  Every  nonstationary  boundary,  be  it  an  agitat¬ 
ing  grid,  propeller  or  moving  walls  of  a  container, 
must  be  heated  to  a  sufficiently  high  temperature  so 
that  frazil  ice  does  not  adhere  to  it  (or  otherwise 
made  ice-phobic).  Otherwise,  frazil  ice  will  accumu¬ 
late  on  the  moving  boundary,  which  can  then  sweep 
it  through  the  fluid  and  lead  to  high  rates  of  ice 
nuclei  formation  and  frazil  ice  dispersion  that  would 
be  artifacts  of  the  particular  apparatus. 

2.  The  primary  fluid  motion  must  persist  for 
several  minutes  without  external  influences.  This 
period  is  required  for  the  growth  of  the  initial  ice 
nuclei  or  crystals,  the  transfer  of  heat  from  them, 
their  rupture,  their  dispersion  and  the  formation  of 
frazil  ice  floes  until  the  sensible  heat  associated  with 


latent  heat  of  fusion  accompanying  ice  growth  is 
absorbed  by  the  supercooling.  This  consideration 
precludes  the  use  of  relatively  short  recirculating 
flumes  for  many  types  of  frazil  ice  experiments,  be¬ 
cause  the  agitation  and  attendant  frazil  ice  forma¬ 
tion  that  may  take  place  in  the  pump  and  during 
flow  passage  through  the  return  pipe  completely 
dominates  the  frazil  ice  formation  that  takes  place 
in  the  working  section  of  the  flume. 

3.  The  apparatus  must  permit  easy  control  of  the 
temperature  or  degree  of  supercooling  of  the  water. 

4.  The  apparatus  must  be  constructed  so  that  a 
controllable  quantity  of  seed  crystals  for  frazil  ice 
nucleation  can  be  easily  introduced  into  the  super¬ 
cooled  water. 

5.  The  water  being  used  in  the  experiments  must 
be  shielded  from  intrusion  of  ice  crystals  from  the 
air  or  other  sources. 

A  turbulence  jar  was  found  to  be  best  suited  to 
the  present  experiments,  which  used  a  turbulence 
jar  similar  to  that  described  by  Mueller  (1978)  and 
Mueller  and  Calkins  (1978).  However,  the  jar  was 
modified  so  that  the  supercooled  volume  of  water 
could  be  isolated  from  the  chilling  elements  that 
cooled  the  water.  The  turbulence  jar  is  shown  in 
Figures  3  and  4.  Figure  3  shows  the  experimental 
setup,  while  Figure  4  is  a  closeup  of  the  turbulence 
jar.  Dimensions  of  the  jar  are  given  in  Figure  5.  The 
metal  sidewall  and  bottom  panels  of  Mueller’s  turbu¬ 
lence  jar,  which  were  used  for  chilling  the  water. 


Figure  4.  Turbulence  jar. 
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Figure  5.  Dimensions  of  the  turbulence  jar. 

frequently  became  covered  by  a  layer  of  boundary- 
fast  ice.  Ice  that  freezes  to  the  jar  boundaries  is  lost 
from  the  frazil  ice  field,  reducing  the  number  of  nu¬ 
cleation  centers  present  in  the  liquid  water  and  there¬ 
by  affecting  the  experimental  results.  For  the  modi¬ 
fied  turbulence  jar,  once  the  desired  temperature  of 
supercooling  was  attained,  the  central  water  com¬ 
partment  was  isolated  from  the  two  chiller  elements. 
The  side  compartments  containing  the  chiller  ele¬ 
ments  were  then  drained  to  further  isolate  the  chiller 
elements  from  the  test  compartment. 

The  water  in  the  turbulence  jar  was  supercooled 
by  circulating  a  chilled  ethylene  glycol  solution 
through  two  chiller  elements.  The  ethylene  glycol 
solution  was  maintained  at  a  temperature  of  -2°C. 

An  oscillating  grid  attached  to  a  scotch  yoke  (see 
Fig.  3)  was  used  to  agitate  the  water  in  the  jar.  The 
grid  was  made  of  6-mm-diameter  tubing  that  was 


Figure  3.  Fxperimental  apparatus. 


bent  into  loops  spaced  30  mm  apart  (see  Fig.  4). 
During  the  phase  of  each  experiment  in  which 
frazil  ice  was  being  formed,  the  gr'd  was  heated  in 
order  to  stop  the  frazil  ice  from  adhering  to  it. 

Early  tests  were  conducted  with  an  unheated 
grid,  but  these  tests  proved  to  be  unsuccessful  be¬ 
cause,  within  a  short  time  after  seeding,  the  frazil 
ice  adhered  to  the  grid.  Thereafter,  as  the  grid  os¬ 
cillated,  it  moved  the  frazil  ice  conglomeration  ad¬ 
hering  to  it  through  the  water  and  in  the  process 
accumulated  more  frazil.  Consequently,  the  mech¬ 
anisms  of  interest -frazil  ice  crystal  dispersion  and 
frazil  ice  floe  and  particle  rupture  by  turbulence- 
were  no  longer  the  dominant  mechanisms  in  the  ex¬ 
periments;  turbulent  diffusion  was  insignificant  com¬ 
pared  to  frazil  advection  by  the  moving  grid.  More¬ 
over,  the  turbulence  produced  by  the  grid  motion 
changed  as  ice  accumulated  on  the  grid  and  altered 
its  effective  dimensions. 

The  rise  of  water  temperature  because  of  heat 
input  from  the  heated  grid  and  from  the  air  surround¬ 
ing  the  turbulence  jar  provided  a  “base-line”  of  tem¬ 
perature  rise  with  time,  as  is  shown  in  Figure  2;  tem¬ 
perature  rise  above  this  line  is  ascribable  to  latent 
heat  release  during  the  formation  of  frazil  ice.  Dur¬ 
ing  the  experiments,  air  temperature  around  the  jar 
was  kept  at  about  1°  to  2°C. 

Much  effort  was  devoted  to  perfecting  the  experi¬ 
ments.  Yet,  although  many  of  the  difficulties  in¬ 
volved  in  conducting  frazil  ice  experiments  were 
overcome,  several  problems  could  not  be  readily 
solved  or  circumvented.  Some  preliminary  experi¬ 
ments  were  also  conducted  in  the  refrigerated  flume 
operated  by  the  I1HR  and  in  a  couette-flow  annulus. 
However,  both  of  these  apparatus  were  found  not 
suitable  for  studying  the  influences  of  water  tempera¬ 
ture  and  turbulence  on  frazil  ice  formation.  The  con¬ 
clusions  derived  from  the  use  of  the  flume  and  the 
couette-flow  annulus  are  briefly  discussed  in  Appen¬ 
dix  A. 

Experimental  procedure 

During  each  experiment,  the  deionized  water  con¬ 
tained  in  the  turbulence  jar  was  supercooled  using 
the  two  chiller  elements  depicted  in  Figure  5.  In 
order  to  prevent  an  illegitimate  nucleation  of  the 
water  by  extraneous  ice  crystals  introduced  from  the 
air  surrounding  the  turbulence  jar,  the  air  tempera¬ 
ture  in  the  ice  room  was  maintained  at  1°  to  2°C. 
Water  temperature  in  the  turbulence  jar  was  registered 
by  means  of  a  fast-response  thermistor  probe.  The 
probe  was  positioned  so  as  to  protrude  horizontally 
about  20  mm  from  one  sidewall  into  the  water  at  its 
mid-depth.  The  probe’s  tip  had  to  extend  at  least  7 
mm  from  the  sidewall  to  be  sufficiently  outside  a 


thin  temperature  boundary  layer  at  the  walls  of  the 
turbulence  jar. 

When  the  water  in  the  turbulence  jar  had  been 
cooled  to  a  temperature  just  below  the  desired  tem¬ 
perature  of  supercooling,  the  central  control  volume 
was  isolated  from  the  chiller  grids  by  the  insertion  of 
the  two  sliding  sidewalls.  At  the  same  time,  the  heat¬ 
ing  element  of  the  agitator  grid  was  activated.  On  at¬ 
taining  the  specified  temperature  of  supercooling,  a 
small  piece  of  ice,  having  a  volume  of  approximately 
0.025  cm3 ,  was  dropped  into  the  water.  After  a  peri¬ 
od  of  about  5  to  10  seconds,  depending  on  the  tem¬ 
perature  of  supercooling  and  the  frequency  of  grid 
oscillation,  frazil  ice  crystals  became  visable.  The 
temperature  of  the  water  and  ice  mixture  was  record¬ 
ed  continuously,  using  a  Gould  strip-chart  recorder, 
until  the  temperature  slightly  exceeded  0°C.  The 
temperatures  of  supercooling  T0  and  frequencies  of 
grid  oscillation  Fg  that  were  used  in  the  study  are 
listed  in  Table  1.  Several  tests  were  conducted  for 
each  combination  of  T0  and  Fg. 

A  series  of  tests  was  conducted  to  measure  the 
growth  rate  of  frazil  ice  particles  with  time  for  vari¬ 
ous  values  of  T0>F%  and  e.  The  frazil  ice  particles 
were  removed  from  the  turbulence  jar  by  means  of  a 
fine-mesh  net  fixed  to  the  end  of  a  spatula.  The 
spatula  was  chilled  to  just  below  0°C  prior  to  its  im¬ 
mersion  into  the  turbulence  jar.  The  spatula  with 
frazil  ice  particles,  on  removal  from  the  jar,  was 
placed  in  a  freezer  box.  A  six-power  optical  compar¬ 
ator  with  a  reticle  graduated  into  0.1 -mm  increments 
was  used  to  measure  the  size  of  the  frazil  ice  particles. 

Table  1 .  Program  of  experiments. 

Initial  supercooling  Grid  oscillation 
temperature  1°  C) _ (cycles/min.) 

-0.05  15.30,60,75,100 

-0.10  15,30,60,75,100 

-0.25  15,30,60,75,100 

-0.55  15,30,45,60,75,100 

-0.60  15,30,45,60,100 

-0.70  15,30,45,60,100 

-0.90 _ 15,30,45,60,75,100 

Values  of  the  turbulence  exchange  coefficient  e 
for  the  agitated  water  in  the  turbulence  jar  were  de¬ 
termined  by  measuring  the  vertical  distribution  of 
suspended  sediment  concentration  in  the  turbulence 
jar.  A  20-mm-thick  layer  of  0.1  -mm-diameter  sand 
was  placed  over  the  base  of  the  jar;  the  action  of  the 
oscillating  grid  caused  some  of  the  sediment  to  be¬ 
come  suspended  throughout  the  depth  of  the  turbu¬ 
lence  jar.  Concentrations  of  suspended  sediment 
were  determined  by  withdrawing  a  volume  of  water 
and  sediment  at  several  elevations  in  the  jar.  A 


narrow-nozzle  pipette,  with  a  capacity  of  10  cm3 , 
was  used  to  extract  the  sample  volumes  of  sediment- 


of  initial  temperature  of  supercooling  and  frequency 
of  grid  oscillation  listed  in  Table  1  are  presented. 
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bearing  water.  The  sediment  contained  in  each 
sample  was  then  weighed.  The  weight  of  sediment 
in  a  sample  with  a  low  sediment  concentration  was 
determined  by  carefully  sluicing  the  sediment  onto 
a  Hat-bottomed  dish  and  then  estimating,  with  the 
aid  of  the  optical  comparator,  its  volume  and  then 
calculating  its  weight.  The  weights  of  sediment  in 
the  samples  with  higher  sediment  concentrations 
were  determined  by  direct  weighing.  The  weight 
concentrations  of  the  samples  were  estimated  from 
the  weights  of  sediments  contained  in  the  sample 
volumes. 


The  temporal  records  of  water  temperature  rise  as 
determined  using  the  analytical  model  are  given  in 
Appendix  C.  The  rise  in  water  temperature  is  propor¬ 
tional  both  to  the  amount  of  frazil  ice  formed  in  the 
turbulence  jar  and  to  the  heat  input  from  the  oscil¬ 
lating  heated  grid.  The  experimental  data  were  used 
to  calibrate  the  analytical  model,  which  was  used  to 
determine  the  rate  of  frazil  ice  formation  as  a  func¬ 
tion  of  water  temperature  and  turbulence. 


Nucleation  of  frazil  ice 

Various  methods  for  seeding  the  supercooled 


RESULTS 

Introduction 

Records  of  the  temporal  rise  of  water  tempera¬ 
ture  for  the  experiments  involving  the  combinations 


water  in  the  turbulence  jar  were  tested.  The  two 
principal  methods  that  were  used  involved  either 
withdrawing  a  specified  volume  (usually  about  20 
cm3 )  using  a  syringe,  exposing  it  to  an  ice  crystal, 
then  injecting  the  fluid  plus  ice  nuclei  back  into  the 
turbulence  jar,  or  dropping  a  “standard”  piece  of  ice 
(about  0.025  cm2  in  volume)  into  the  turbulence  jar. 


For  each  method  there  was  no  difference  in  the 
manner  of  frazil  ice  formation,  provided  that  the 
water  and  ice  in  the  syringe  were  not  vigorously 
agitated;  when  this  happened,  a  larger  amount  of 
frazil  ice  nuclei  could  be  introduced  into  the  turbu¬ 
lence  jar  and  thus  lead  to  a  slightly  greater  rate  of 
frazil  ice  formation.  For  most  of  the  experiments, 
the  turbulence  jar  was  initially  seeded  by  dropping 
in  a  standard  piece  of  ice. 

As  in  the  experimental  studies  by  Mueller  (1978) 
and  others,  as  well  as  in  the  field  observations  re¬ 
ported  by  Osterkamp  (1978b),  this  study  found  that 
there  was  no  spontaneous  or  homogeneous  nuclea- 
tion  of  the  supercooled  water.  In  all  experiments, 
seeding  was  required  to  initiate  frazil  ice  formation. 
The  formation  of  frazil  ice  in  the  turbulence  jar  is 
illustrated  in  the  series  of  photographs  presented  in 
Figure  6.  Figure  6a  depicts  the  water  being  cooled 


to  a  temperature  of  -0.55  C  at  a  frequency  of  grid 
oscillation  of  75  cycles/min.  Small  disks  of  frazil 
ice  are  forming  in  the  supercooled  water,  20  seconds 
after  nucleation  (Fig.  6b).  Some  50  to  70  seconds 
after  nucleation  (Fig.  6c  and  d  respectively),  the 
frazil  ice  particles  have  increased  in  size  to  a  mean 
diameter  of  about  6  to  7  mm. 

Subsequent  to  the  seeding  of  the  water,  a  brief 
period,  typically  5  to  10  seconds,  elapsed  before 
frazil  ice  particles  appeared  in  the  turbulence  jar. 
Once  a  few  frazil  jce  ice  particles  had  appeared,  the 
particles  multiplied  rapidly.  The  rapid  increase  in 
the  number  of  frazil  ice  particles  formed  can  be  as¬ 
cribed  to  the  process  of  collision  breeding;  frazil  ice 
particles,  by  colliding  with  the  solid  boundaries  of 
the  turbulence  jar,  the  oscillating  grid  or  with  each 
other,  produce  an  increasing  number  of  ice  nuclei. 


-  t 
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Influences  of  turbulence  on  frazil  ice  formation 

As  the  turbulence  intensity  of  a  flow  increases, 
the  rate  of  frazil  ice  formation  increases.  The  in¬ 
crease  in  the  rate  of  frazil  ice  formation  is  primari¬ 
ly  caused  by  both  the  increased  rate  of  heat  trans¬ 
fer  between  the  frazil  ice  particles  and  the  surround¬ 
ing  water  and  an  increased  rate  of  floe  rupture, 
which  provides  more  ice  nuclei.  As  shown  by  the 
experimental  data,  the  influence  of  turbulence  on 
the  rate  of  frazil  ice  formation  appears  to  be  more 
pronounced  for  lower  temperatures  of  supercooling. 
This  result  can  be  observed  from  the  records  of  ex¬ 
perimental  data  that  are  presented  in  Figure  7. 


The  solid  line  curves  in  Figure  7  were  fitted  to 
temperature  measurements  that  were  recorded  for 
several  experiments  at  constant  values  of  T0  and  fre¬ 
quency  of  grid  oscillation.  The  dash-dot  curves  are 
plotted  in  order  to  indicate  the  range  of  the  experi¬ 
mental  data.  Although  the  trends  for  temperature 
rise-that  is,  frazil  ice  formation-indicate  that  the 
rate  of  frazil  ice  formation  generally  increases  with 
greater  values  of  grid  oscillation  frequency,  a  clear 
trend  is  not  always  apparent,  especially  for  values  of 
Ta  equal  to  -0.1 0°C  and  -0.05°C. 

It  would  appear  from  the  experimental  data  that 
the  rate  of  frazil  ice  formation  may  not  be  quite  as 


TIME  SINCE  SEEDING  (SECONDS) 

0  20  40  60  80  100  120  140 


a.  Ta  =  -0.90°  C. 


» 


b.  T o=-0.70°C. 


Figure  7.  Experimental  data.  The  solid  curves  are  plotted  for  the  mean  results  for  each  grid  frequency,  and  the 
dash-dot  lines  indicate  the  range  of  the  data  (CPM -cycles  per  minute ). 
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Figure  7  (cant ’d).  Experimental  data.  The  solid  curves  are  plotted  for  the  mean  results  for  each  grid  frequency, 
and  the  dash-dot  lines  indicate  the  range  of  the  data  (CPM -cycles  per  minute). 
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Figure  7  (cont’d). 


sensitive  to  turbulence  intensity  as  has  been  suggested 
in  recent  literature  (Mueller  1978,  Mueller  and  Calkins 
1978,  Martin  1981).  Rather,  turbulence  intensity  may 
indirectly  exert  a  greater  influence  on  the  rate  and 
quantity  of  frazil  ice  formation  by  influencing  the 
temperature  to  which  a  body  of  water  can  be  super¬ 
cooled.  When  Ta  is  equal  to  either  -0.1 0°C  or 
-0.05°C,  the  relationship  between  the  rate  of  tem¬ 
perature  rise  and  frequency  of  grid  oscillation  ap¬ 
pears  to  weaken.  In  both  Figures  7f  and  g,  bands 
enclosing  the  experimental  data  are  shown  and  the 
overlap  of  some  of  the  temperature-time  curves 
is  indicated.  Although  the  rate  of  frazil  ice  forma¬ 
tion  increases  slightly  with  increased  level  of  turbu¬ 
lence,  the  trend  is  not  sufficiently  marked  that  the 
vagaries  of  experimentation  cannot  mask  it. 

The  rates  of  frazil  ice  formation  when  Ta  = 

-0.55°C  are  shown  in  Figure  7d  for  two  temperature 
settings  of  the  heated  grid.  Two  settings  for  grid 
heating  were  used  in  order  to  test  its  effect  on  the 
rate  of  frazil  ice  formation.  As  is  subsequently  in¬ 
dicated  in  Figure  10,  the  heat  input  into  the  jar 
from  the  grid  did  not  have  a  pronounced  effect  on 
the  rate  of  frazil  ice  formation  for  larger  tempera¬ 


tures  of  supercooling.  The  heat  input  was,  however, 
reduced  for  values  of  T0  =  -0.25 °C,  -0. IO°C  and 
-0.05°C,  with  the  intention  that  its  contribution  to 
the  heating  of  the  jar  would  not  dominate  that  of 
the  latent  heat  of  fusion  resulting  from  frazil  ice 
growth.  Nevertheless,  it  was  found  (using  the  analyt¬ 
ical  model)  that  when  Ta  =  -0.1 0°C  and  -0.05°C  the 
heat  input  from  the  grid  affected  the  experimental 
data.  Because  the  heat  input  T  from  the  grid  could 
be  accounted  for  in  the  analytical  model,  eq  17 
through  19,  the  experimental  data  could,  however, 
be  used  to  both  validate  and  calibrate  the  model. 

From  the  comparison  of  experimental  data  and 
calculated  results  obtained  using  eq  17  through  19, 
and  using  equivalent  values  of  T ,  it  was  found  that 
during  an  experiment  the  heat  transfer  coefficient 
a  varies  from  an  initial  value  of  a0  to  approximate¬ 
ly  one-half  of  its  value  (0.5  aG).  Correlation  be¬ 
tween  experimental  and  analytical  data  led  to  selec¬ 
tion  of  the  values  of  the  constants  given  in  eq  19 
with  the  following  results: 


a  =  0.5  a0  (I  +  e  <4°  c,>). 


The  decreasing  value  of  the  heat  transfer  coeffi¬ 
cient  a  with  time  is  attributable  both  to  the  increas¬ 
ing  viscosity  of  the  water-ice  mixture  as  the  water 
temperature  in  the  turbulence  jar  rises  and  to  the 
interference  of  the  frazil  ice  particles  with  the  heat 
transfer  process.  Frazil  ice  particles  suspended  in 
the  fluid  act  to  damp  the  turbulence  eddies  in  the 
flow,  thus  leading  to  a  reduction  in  turbulence  in¬ 
tensity  and  resulting  in  a  decline  in  the  value  of  the 
heat  transfer  coefficient.  Another  factor  contribut¬ 
ing  to  the  reduction  in  the  rate  of  heat  transfer  from 
the  frazil  ice  particles  is  that  the  particles  agglutinate 
or  sinter  (Martin  1981 )  and  agglomerate  as  floes, 
thereby  reducing  the  surface  area  of  ice  that  is  ex¬ 
posed  to  the  water. 

Predicted  records  of  temperature  rise  with  time 


are  presented  in  Appendix  C.  These  records  were 
obtained  from  the  numerical  integration  of  eq  17 
and  18  with  a  expressed  as  eq  28.  In  Figure  Cl  the 
same  values  of  T  were  used  as  were  used  in  obtain¬ 
ing  the  experimental  data,  but  in  Figure  C2,  T  was 
set  equal  to  zero.  Comparison  of  Figures  Cl  and  C2 
for  t  =  -0.1 0°C  and  -0.05°C  shows  the  degree  to 
which  grid  heating  marked  the  experimental  data  for 
these  values  of  T0. 

The  frequency  of  grid  oscillation  was  related  to 
turbulence  exchange  coefficient  e  in  the  following 
way.  The  variation  with  depth  of  the  concentration 
of  suspended  fine  sand  particles  in  the  turbulence 
jar  for  several  frequencies  of  grid  oscillation  was 
plotted  as  is  shown  in  Figure  8.  From  Figure  8, 
Figure  9  was  plotted  to  obtain  (in  the  manner 
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Figure  8.  Distributions  of  suspended  sediment  concentration  with  depth  in  the  turbulence 
jar  versus  frequency  of  grid  oscillation  Fg  (CPM -cycles  per  minute). 
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Figure  9.  Relationship  between  frequency  of  grid  oscillation  Fy  and  turbu¬ 
lence  exchange  coefficient  e  (CPM -cycles  per  minute). 
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described  in  the  Elements  of  Turbulence  section) 
the  relationship  between  frequency  of  grid  oscilla¬ 
tion  Eg  and  the  corresponding  value  of  turbulence 
exchange  coefficient  e.  It  should  be  noted  that  the 
values  of  e  obtained  for  the  jar  are  gross  values.  Be¬ 
cause  of  the  prismatic  form  of  the  turbulence  jar 
and  the  shape  of  the  oscillating  grid,  values  of  e 
would  vary  somewhat  from  location  to  location 
within  the  jar. 

The  average  temporal  rate  of  temperature  rise 
because  of  frazil  ice  formation  in  the  tuibulence  jar 


B  razil’  multiplied  by  the  specific  heat  (or  heat  capa¬ 
city )  of  water  at  0°C,  SH  =4.218  kJ/kg°C,  is  plotted 
in  Figure  10  against  the  grid  oscillation  Fg.  The  or¬ 
dinate  in  Figure  10  is  equivalent  to  the  rate  of  heat 
release  per  unit  weight  of  water  because  of  frazil  ice 
formation,  that  is 

Qw  -  p  l  -S'h  B  razil 
or 

Gw 

-pV  -  S»  TVrazi.  (29) 
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Figure  10.  Average  temporal  rates  of  frazil  ice  formation  versus  fre¬ 
quency  of  grid  oscillation  F  initial  heat  transfer  coefficient  a(),  and 
turbulence  exchange  coefficient  e.  The  curves  are  plotted  through 
the  analytically  determined  data  (CPM  cycles  per  minute). 


where  =  average  temporal  rate  of  temperature 
rise  because  of  frazil  ice  formation.  The  reader 
should  bear  in  mind  that  the  effect  of  the  tempera¬ 
ture  rise  from  the  heated  grid  has  been  removed 
from  Figure  10.  Plotted  as  an  alternate  abscissa  in 
Figure  10  are  values  of  the  turbulence  exchange  co¬ 
efficient  e  determined  for  the  frequencies  of  grid  os¬ 
cillation  F  .  Also  plotted  in  Figure  10  are  the  values 
of  SH  4azi.  as  predicted  using  the  analytical  model 
of  the  frazil  ice  formation.  The  predicted  values  of 
•SH  ^Frazil  are  pl°tted  against  the  initial  value  of  the 
heat  transfer  coefficient  aG.  The  abscissa  of  Figure 
10  is  composed  of  three  scales,  F  ,  e  and  ao.  The 
a0  scale  was  set  so  that  the  experimental  and  the 
analytical  data  coincided  and  the  analytical  model 
was  thereby  calibrated.  The  curves  shown  in  Figure 
10  are  drawn  through  the  values  of  SH  /^.-razi,  as  de¬ 
termined  using  eq  17,  18  and  28.  For  the  same 
values  of  F  ,  or  equivalent  a0,  both  the  analytically 
predicted  and  the  experimentally  determined  values 
of  ^Frazil  are  g*ven-  Although  the  experimental 
data  exhibit  some  scatter,  the  analytically  predicted 
data  overlay  the  experimental  data  fairly  well.  Both 
the  experimental  and  the  calculated  data  are  associ¬ 
ated  with  equivalent  values  of  grid  heating  f . 

It  is  apparent  from  Figure  10  that  values  of  tur¬ 
bulence  exchange  coefficient  e  ranging  from  0  to 
about  20  cm2/s  correspond  to  initial  values  of  heat 
exchange  coefficient  aQ  ranging  from  1  to  14  W/°C 
m2.  Figure  10  indicates  that  the  analytical  model 
adequately  embodies  the  physics  of  frazil  ice  forma¬ 
tion;  the  experimental  and  analytical  data  coincide 
quite  well. 

On  the  basis  of  the  results  given  in  Figure  10,  it 


can  be  concluded  that  the  average  rate  of  frazil  ice 
formation  increases  for  higher  levels  of  turbulence 
in  a  flow.  However,  the  influence  of  turbulence  de¬ 
creases  with  increasing  values  of  initial  temperature 
of  supercooled  water;  for  values  of  T0  less  than 
about  0.10°C,  the  level  of  turbulence  does  not  ap¬ 
pear  to  have  a  major  effect  on  the  rate  of  ice  forma¬ 
tion  in  a  parcel  of  supercooled  water  that  is  not 
bounded  by  a  free  surface. 

Water  temperature 

It  is  evident  from  both  the  experimental  and  ana¬ 
lytical  results  that  frazil  ice  forms  more  rapidly,  and 
in  greater  amounts,  in  colder  supercooled  water. 
These  results  are  graphically  illustrated  in  Figure  1 1 , 
in  which  the  increase  in  water  temperature  attrib¬ 
utable  to  the  formation  of  frazil  ice  is  compared  for 
two  initial  temperatures  of  supercooling,  T0  = 
-0.55°C  and  -0.90°C.  For  the  same  frequency  of 
grid  oscillation,  or  level  of  turbulence,  the  initially 
colder  water  is  warmed  at  a  greater  rate  by  the 
greater  amount  of  frazil  ice  that  forms.  The  afore¬ 
mentioned  variation  of  the  amount  of  frazil  ice 
formation  is  also  depicted  in  Figure  12,  in  which 
the  amount  of  frazil  ice  formed  as  estimated  using 
the  analytical  model-that  is  using  eq  17,  18  and  28- 
is  plotted  against  time  for  parcels  of  water  with  the 
same  value  of  turbulence  exchange  coefficient  e. 

An  inspection  of  eq  17  reveals  that  the  rate  of 
temperature  rise  (i.e.  the  rate  of  frazil  ice  formation) 
is  a  quadratic  in  T.  Consequently,  the  lower  the 
temperature  of  supercooling  T0,  and  hence  the  ini¬ 
tial  values  of  T,  the  greater  the  rate  of  frazil  ice 
formation. 
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Figure  11.  Comparison  of  rates  of  temperature  rise  attributable  to  frazil  ice 
formation  for  two  initial  temperatures  of  supercooled  water  and  similar 
levels  of  turbulence  (CPM -cycles  per  minute). 
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Figure  12.  Comparison  of  rates  of  frazil  ice  formation  for  tsvo  initial  tempera¬ 
tures  of  supercooled  water  and  similar  levels  of  turbulence. 


Influences  of  water  temperature  and  turbulence 
on  the  concentration  of  frazil  ice 

The  concentration  of  frazil  ice  as  a  function  of 
time  for  several  temperatures  of  supercooling  and 
initial  values  of  the  heat  transfer  coefficient  are 
given  in  Figure  13.  The  curves  in  Figure  13  were 
uetermined  using  the  analytical  model,  that  is,  eq 


1 7  with  f  set  as  zero  and  eq  1 8  and  28  calibrated 
to  the  experimental  results  (Fig.  10).  The  average 
temporal  rates  of  frazil  ice  formation,  calculated 
as  the  final  volume  of  frazil  ice  divided  by  the  time 
for  the  water  temperature  to  reach  0°C,  are  plotted 
with  aQ  and  e  and  are  presented  in  Figure  14. 


a.  T0  =  -0.90° C. 


Figure  13.  Temporal  variation  in  frazil  ice  concentration  for  dif¬ 
ferent  values  of  initial  heat  transfer  coefficient,  a0  (T  =  0). 
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Figure  14.  Variation  of  the  average  rate  of  frazil  ice  formation  with  tempera¬ 
ture  of  supercooling  T0,  initial  value  of  heat  transfer  coefficient  a0  and  tur¬ 
bulence  exchange  coefficient  e. 


With  larger  initial  temperatures  of  supercooling 
Ta ,  both  the  rate  of  frazil  ice  formation  and  the 
final  concentration  of  frazil  ice  increase.  The  tem¬ 
poral  rate  of  frazil  ice  formation  (measured  as  volu¬ 
metric  concentration)  also  increases  with  increasing 
values  of  turbulence  exchange  coefficient,  that  is, 
with  increasing  turbulenc  5  intensity.  It  is  evident 
from  Figure  14  that  the  following  relationship  exists 
between  c,  e  and  T0 : 

c  =  e"  f(Ta)  (30) 

where  n  ^  1 .0  for  all  values  of  T0.  In  other  words, 
the  average  rate  of  frazil  ice  formation  is  directly 
proportional  to  the  turbulence  intensity  of  the 
fluid  in  which  it  forms.  The  average  rate  of  frazil 


ice  formation  is  strongly  related  to  the  initial  tem¬ 
perature  of  supercooling  T0. 

No  measurements  were  made  of  the  actual  con¬ 
centration  of  frazil  ice  that  formed  in  the  turbulence 
jar.  A  review  of  the  present  techniques  for  measur¬ 
ing  frazil  ice  concentration  (for  example,  that  pro¬ 
posed  by  Tsang  1977)  showed  that  these  techniques 
require  further  development  before  they  can  be  used. 
The  concentration  of  frazil  ice,  for  this  study,  was 
determined  analytically  by  way  of  the  heat  produced 
and  the  attendent  water  temperature  rise  accompany¬ 
ing  frazil  ice  formation.  The  analytical  model  indi¬ 
cates  that  the  final  concentration  of  frazil  ice  for  a 
given  temperature  of  supercooling  T0  is  independent 
of  the  turbulence  regime  of  the  flow.  Turbulence 
affects  the  rate  of  frazil  ice  formation  but  not  its 
final  concentration.  Turbulence  may  indirectly 
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TEMPERATURE  OF  SUPERCOOLING  (-°C) 

Figure  15.  Variation  of  the  final  concentration  of  frazil  ice 
with  temperature  of  supercooling  T0. 


influence  the  amount  of  frazil  ice  formed  by  enabl¬ 
ing  a  body  of  water  to  become  supercooled  to  a 
lower  temperature  prior  to  nucleation.  The  final 
values  of  frazil  ice  concentration  cmax  are  plotted 
in  Figure  15  against  Ta.  In  keeping  with  the  formu¬ 
lation  of  the  analytical  model,  the  values  of  cmax 
are  directly  proportional  to  T0.  On  the  basis  of  the 
analytical  model,  calibrated  with  the  experimental 
data,  the  following  relationship  can  be  proposed  for 


cmax  =  “0-014  T0 


where  cmax  =  maximum  volumetric  concentration 
of  frazil  ice  produced  in  supercooled  water  and  Ta  = 
minimum  temperature  below  freezing  (°C)  of  the 
supercooled  water.  Equation  3 1  is  essentially  the 
theoretical  limit  for  the  concentration  of  ice  formed 
in  a  parcel  of  supercooled  water.  This  limit  can  be 
estimated  by  equating  the  quantity  of  ice  that  is  re¬ 
quired  to  raise  the  water  temperature  to  0°C  by  way 
of  release  of  latent  heat  of  fusion  Lx.  That  is, 


=  -0.014  T 


where  p  and  p(  =  the  densities  of  liquid  water  and 
ice  at  the  temperature  of  supercooling  T0,  SH  = 


specific  heat  of  water,  4.218  kJ/kg°C  and  Li  =  latent 
heat  of  fusion  of  ice,  333.6  kJ/kg.  For  river  water 
supercooled  to  about  -0.1°C,  eq  31  gives  a  maximum 
volumetric  concentration  of  frazil  ice  of  0.14%.  This 
value  is  approximately  one  quarter  of  the  value  that 
was  reported  by  Tsang  (1982),  who  cites  a  maximum 
volumetric  concentration  of  0.55  %  for  rivers  (from 
Ontario  Hydro  1970).  The  disparity  in  the  two  values 
may  have  one  of  several  causes.  The  field  value,  cited 
by  Tsang,  was  evaluated  from  samples  taken  from  the 
upper  portion  of  flow  in  the  Niagara  River.  Observed 
concentrations  of  frazil  ice  are  likely  to  be  higher  in 
the  upper  portion  of  a  flow  because  of  the  influence  of 
buoyancy  that  causes  frazil  ice  particles  to  rise  toward 
the  surface.  An  additional  possible  cause  is  that  the 
river  may  continuously  lose  heat  to  frigid  air,  thereby 
enabling  frazil  ice  formation  to  continue.  Larger  con¬ 
centrations  of  frazil  ice  may  occur  when  the  particles 
agglutinate  and  locally  agglomerate  into  slush  ice; 
slush  ice,  however,  forms  as  a  consequence  of  the 
transport  and  mechanical  concentration  of  frazil  ice 
by  the  flow  and  is  only  indirectly  dependent  on  the 
formation  of  frazil  ice  from  a  supercooled  parcel  of 
water. 

Frazil  ice  particle  shape  and  size 

Measurements  of  the  particle  size  of  frazil  ice 
showed  that  maximum  particle  size,  both  in  terms 
of  the  length  of  the  longest  axis  and  the  surface  area, 


increases  with  decreasing  initial  water  temperature  the  more  rapid  rate  of  frazil  ice  formation  and  long- 

Ta  and  decreases  with  increased  turbulence  intensity.  er  period  for  growth.  At  the  higher  turbulence  inten 

Frazil  ice  particles  as  single  crystals  and  as  agglu-  sities,  a  limiting  size  of  frazil  ice  particles  is  reached 

tinated  or  sintered  laminations  of  crystals  were  ob-  because  the  turbulence  eddies  act  to  buffet,  flex  and 

served  and  measured.  Size  distributions  of  frazil  ice  break  those  particles  that  get  large.  Measured  values 

particles  were  not  determined.  Larger  frazil  ice  par-  of  maximum  particle  length  or  diameter  (as  defined 

tides  can  form  in  initially  colder  water  because  of  in  the  sketches  in  Fig.  16)  of  observed  frazil  ice 
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Figure  16.  Summary  of  sizes  and  shapes  of  frazil  ice  particles  observed  when 
temperature  of  supercooling  was  -0.35° C. 


Figure  1 7.  Measured  values  of  maximum  diameter  of  frazil  ice  particles  versus  water  temperature 
normalized  by  temperature  of  supercooling  T„  (CPM  -cycles  per  minute). 
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Figure  18.  Final  maximum  length,  or  diameter,  of 
frazil  ice  particles  versus  temperature  of  supercool¬ 
ing  T  (CPM  cycles  per  minute). 


particles  are  plotted  in  Figure  1 7  against  water  tem¬ 
perature  that  is  normalized  by  initial  temperature 
of  supercooling  T/T0.  The  final  maximum  lengths 
of  frazil  ice  particles  are  plotted  in  Figure  18  against 
temperature  of  supercooling,  with  turbulence  ex¬ 
change  coefficient  e  as  a  third  parameter.  The  trend 
of  increasing  frazil  ice  particle  size  with  decreasing 
level  of  turbulence  is  evident  in  Figure  18. 

When  the  temperature  of  supercooling  was  less 
than  about  -0.1 0°C,  the  frazil  ice  particles  had  the 
form  of  disks  with  maximum  diameter  of  approxi¬ 
mately  2  to  3  mm  and  thickness  of  about  0.05  mm. 
For  water  supercooled  to  initially  lower  tempera¬ 
tures,  the  frazil  ice  particles  had  sufficient  time  to 
grow  to  a  larger  size  and  to  become  more  irregular 
in  shape  (see  Fig.  16).  Occasionally,  the  larger 
frazil  ice  pai  tides  would  have  smaller  particles 
(single  crystals)  about  1  mm  in  diameter  fused  or 
sintered  to  their  surfaces. 


CONCLUSIONS 


A  study  involving  experimentation  coupled  with 
an  analytical  model  of  frazil  ice  formation  in  a  super¬ 
cooled  parcel  of  water  was  conducted  in  order  to  in¬ 
vestigate  the  influences  of  turbulence  and  water  tem¬ 
perature  on  frazil  ice  formation.  The  following  con¬ 
clusions  were  d~awn  from  the  study. 

Ice  fragments  or  crystals  are  required  to  seed 
water  that  has  been  supercooled  to  temperatures 
as  low  as  -0.9°C.  Chilled,  nonaqueous  particles 
placed  in  supercooled  water  do  not  cause  frazil 
ice  to  develop. 

The  most  important  determinant  of  frazil  ice  for¬ 
mation  is  the  initial  temperature  of  supercooling. 
Both  the  formation  rate  and  the  concentration  of 
frazil  ice  increase  with  lower  temperatures  of  super¬ 
cooling.  A  simple  relationship  exists  between  the 
temperature  of  supercooling  and  the  maximum  con¬ 
centration  of  frazil  ice  that  can  be  formed  in  a  vol¬ 
ume  of  supercooled  water. 

The  rate  of  frazil  ice  formation  is  directly  propor¬ 
tional  to  intensity  of  turbulence  of  the  fluid  in 
which  it  is  forming.  The  turbulence  characteristics 
of  a  flow  affect  the  rate  of  frazil  ice  formation  by 
governing  the  temperature  to  which  the  flow  can  be 
initially  supercooled,  by  influencing  the  processes 
of  heat  transfer  from  the  crystals  to  the  surrounding 
water,  and  by  promoting  collision  nucleation  and 
floe  rupture  and  thereby  increasing  the  number  of 
nucleation  sites. 

Results  from  the  analytical  model  suggest  that 
the  final  concentration  of  frazil  ice,  for  a  given  tem¬ 
perature  of  supercooling,  is  independent  of  the  tur¬ 
bulence  intensity  of  flow. 

Larger  frazil  ice  particles  (crystals  and  platelets 
of  fused  crystals)  form  in  water  supercooled  to  low¬ 
er  temperatures,  but  particle  size  decreases  with  in¬ 
creasing  turbulence  intensity.  Large  frazil  ice  parti¬ 
cles,  having  a  major  platelet  diameter  of  about  1 0 
mm  and  a  thickness  of  about  0.1  mm  or  less,  tend  to 
be  broken  when  jostled  by  flow  eddies. 

For  temperatures  of  supercooling  of -0.1 0°C  and 
less,  as  are  typically  recorded  for  supercooled  water 
in  rivers  and  streams,  the  maximum  length  of  the 
frazil  ice  disks  was  found  to  be  about  2  mm. 
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APPENDIX  A:  PRELIMINARY  FRAZIL 
ICE  EXPERIMENTS 

Flume  experiments 

A  series  of  experiments  was  conducted  using  the 
Iowa  Institute  of  Hydraulic  Research  (1IHR)  40-ft- 
long  ( 1 2  m)  refrigerated  flume  that  used  several 
means  of  ice  crystal  seeding  to  start  nucleation.  The 
object  was  to  obtain  major  growth  of  frazil  ice  dur¬ 
ing  the  passage  of  the  flow  through  the  labyrinth 
channel  formed  by  placing  two  partitions  in  the  2-ft- 
wide  (0.6  m)  flume  to  form  a  120-ft-long  (37  nt), 
8-in.-wide  (0.2  m)  channel. 

In  the  earlier  experiments  with  seeding  near  the 
upstream  end  of  the  channel,  it  was  repeatedly  ob¬ 
served  that  little  frazil  was  formed  along  the  120-ft- 
long  (37  m)  channel,  but  frazil  ice  was  formed  and 
the  supercooling  exhausted  after  the  flow  then 
passed  through  the  pump  and  the  return  pipe.  There¬ 
after,  various  attempts  were  made  to  nucleate  the 
supercooled  flow  near  the  downstream  end  of  the 
channel.  After  the  passage  of  the  seed  ice  through 
the  pump  and  the  return  pipe,  a  large  amount  of 
frazil  ice  appeared  at  the  upstream  end  of  the  flume. 
Apparently,  the  agitation  of  passage  through  the 
pump  and  return  pipe  led  to  formation  of  sufficient 
frazil  to  exhaust  the  supercooling  before  the  flow 
entered  the  open-channel  portion  of  the  flow  circuit. 

In  an  attempt  to  obtain  a  more  rapid  formation 
of  frazil  ice,  experiments  weie  made  in  which  the 
supercooled  flow  was  seeded  with  a  “slug”  of  frazil 
ice,  with  the  objective  of  providing  more  nucleation 
centers.  This  was  done  by  supercooling  water  in  the 
turbulence  jar  (described  in  the  Analytical  Model 
section),  nucleating  it  with  ice  and  injecting  a  por¬ 
tion  of  the  resultant  frazil-water  mixture  into  the 
flow  at  the  upstream  end  of  the  channel.  During 
repeated  runs  there  was,  as  frazil  passed  a  tempera¬ 
ture  probe,  a  temporary  rise  in  local  temperature  in 
the  downstream  direction  along  the  channel.  After 
careful  evaluation  of  this  seemingly  encouraging  re¬ 
sult,  it  was  suspected  that  much  of  the  temperature 
rise  resulted  from  diffusion  of  the  0°C  liquid  water 
introduced  with  the  frazil  ice  near  the  upstream  end 
of  the  flume.  Fortunately  (or  unfortunately)  this 
suspicion  was  confirmed  when  only  frazil  ice  parti¬ 
cles  from  the  turbulence  jar,  instead  of  the  frazil- 
liquid  water  mixture,  was  used  for  nucleation.  In 


these  experiments  there  was  no  appreciable  amount 
of  frazil  ice  formation  or  temperature  rise  along  the 
channel.  The  first  sizable  quantity  of  frazil  that  ap¬ 
peared  was  observed  at  the  upstream  end  of  the 
channel,  after  passage  of  the  seeded  flow  through 
the  pump  and  the  return  pipe. 

So,  the  partitions  were  removed  from  the  flume 
and  experiments  were  made  in  a  high  velocity, 
shallow  supercritical  flow  in  an  effort  to  enhance 
the  rates  of  floe  formation  and  loss  of  supercooling. 
However,  any  increase  in  these  effects  was  more  than 
offset  by  the  reduced  residence  time  of  the  nucleated 
water  in  the  working  section  of  the  flume. 

In  summary,  the  IIHR  flume  proved  to  be  un¬ 
suitable  for  conducting  experiments  to  accomplish 
the  goals  of  the  study. 

Couette-flow 

In  the  quest  for  an  apparatus  that  would  give  an 
endless  turbulent  shear-layer  flow,  which  is  virtually 
free  of  the  end  effects  that  plagued  the  refrigerated 
flume,  a  series  of  experiments  was  conducted  using 
a  couette-flow  annulus.  The  annulus  was  made  from 
0.5-in. -diameter  (1 .3  cm)  tubing  and  had  a  major 
diameter  of  about  8  in.  (20  cm).  It  was  equipped 
with  two  short  radial  nipples  through  which  the  an¬ 
nulus  could  be  filled  with  water  and  nucleated.  The 
annulus  was  attached  to  a  horizontal  axle  through 
its  principal  center  of  curvature  and  the  axle,  in  turn, 
was  attached  to  a  variable-speed  motor.  The  lower 
one-third  of  the  annulus  was  immersed  in  a  tempera¬ 
ture-controlled  salt  water  ice  bath.  The  annulus  was 
filled  to  mid-depth  with  supercooled  water,  then 
nucleated  with  ice  crystals  through  one  of  the  nip¬ 
ples. 

The  frazil  ice  that  formed  subsequently  was 
found  to  adhere  tenaciously  to  the  walls  of  the 
annulus  and  was  then  swept  through  the  sheared 
flow  in  the  lower  part  of  the  annulus  as  it  rotated. 
After  a  few  experiments  it  was  concluded  that  this 
“forced  advection”  of  the  frazil  ice  through  the 
liquid  was  having  the  same  effect  on  frazil  forma¬ 
tion  as  the  transport  of  the  frazil  through  liquid 
water  in  the  early  turbulence  jar  experiments  with 
an  unheated  grid,  and  was  therefore  completely 
dominating  the  normal  particle  diffusion  and  floe 
growth  process.  The  couette-flow  annulus  experi¬ 
ments  were  abandoned. 
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APPENDIX  B:  LISTING  OF  COMPUTER  PROGRAM  FOR 
CALCULATION  OF  FRAZIL  ICE  FORMATION 


//  EXEC  EORTIICLG,REGION=20  0K 
//FOR  r  .  SYS  IN  DD  * 

C  ALPI  l=!  !EAT  TRANSFER  COEFFICIENT 

C  RO=DENSITY  OF  WATER 

C  R 01 “DENSITY  OF  ICE 

C  A I “LATENT  HEAT  OF  FUSION  OF  ICE 

C  DI-CHAR  LENGTH  OF  ICE  PARTICLE ( THICKNESS > 

C  S!’=Sr'ECIEIC  HEAT  OF  ICE 

C  TR -  R ATE  OF  HEATING  BY  GRID 

C  TO=INIT IAL  TEMPERATURE  OE  SUPERCOOLED  WATER 

C  TM-EQUILIDRIUM  TEMPERATURE  OF  ICE  AND  WATER 

C  TEMP -WATER  TEMPERATURE  AT  ANY  TIME 

C  T=T I ME, IN  SECONDS 

C  9=9 GLUME  OF  TUBULENCE  JAR 

C  J  IS  A  COUNTER 

DIMENSION  ALP II  (  260  )  , CO ( 260  )  ,  AD  ( 260  )  ,  TC  (  21.0  >  ,  B  <  260  )  ,  It 
DIMENSION  Cl (260 ) , T ( 260 ) , X ( 260 ) , TEMP ( 260 ) , A ( 260 )  ,  At  <2 
DOODLE  PRECISION  T , f , X , TEMP , A , E , Q , AD , D , ALPI I , TC , CO 
COMMON  RQI ,  RO ,  A.I ,  BI , TR , TO , SI  I , V , D , TM 

C 

N 1 =NQ .  or  ALPHA'S 
N2=N0.  OF  TO'S 
RE AD (5,. 3)  N1,N2 
R E AD ( S , 4 )  ( Ai  ( !<  1 )  ,  K  1  =  1 , Ni  > 

R EAD  (0,4)  ( A2  ( K 2 )  ,  l<2=  i  ,  N2 ) 

FORMAT ( 415 ) 

4  FORMAT (6F10 .2) 

DO  2000  IJ-1,2 
DO  1000  1 1  =::  1  ,  N2 
DO  500  II=i,Ni 

C  INITIALISE  T , TEMP ,  X  >  A  ,  M 

N=i 
J  =  1 
M  =  0 

T0“A2 (II) 

TM- 0 . 00 
1 1=2 . 0 

II"  <  II .  GT  .  5)  11=0. 0 
G 

T ( J  >  =  0 .000 
TEMP ( J ) =TO 
X  (  J  >  =  TEMP  ( J  )  -TM 
A ( J )  =  0 .00  02 
y=0  004128 
LO(..r>  =  0  .  0 0 
c 

C  COMPUTE  STEP  CONSTANTS  IT,  C,  I) , 

ALPI  I  ( J) =Ai(II> 

ir(IJ.EQ.i)  WR ITE ( 6 , 5 )  TO,ALPII(l) 

R 0=997 . 0 
R OI =9 17. 00 
A I =333  60 
SIT =4 . 2177 
DI=0 .0001 
T R = 0  00020 

If (  II  EQ  1 )  I R  = . 00042 
ir(Il.E0.2)  TR=  00085 


w  V  »  K 
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EQ  5) 

IR®  00071 
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1  K  U"l|(,2'> 

I  r 

ai 

:  u  s> 

;  y.  (;  o  ii  i  ■' '• 

•  V 

1  1 

( 1 1 

IQ  6) 

TR-  ooo 

'J. 

IT 

■:  1 1 

E  Q  7) 

1  l\  —  1U  0  46 

A 

D( 

J  1  -A 

LIT  0  J  > 

X  I  R  /  (  H  O I X  A  l  X  V 

I  ) 

1 

C  JT  • 

U '  J  )  /  1  R 

D  ; 

IR 

M 

<  J )  = 

:  ALIM  1  ( 

J  )/  (ROXGI! 

#V  >  ' 

X  ( A  < 

J  >  iROX'TXV 

X ( TM  TO)/ 

(  ROT  X  A  I  x.P T  ) 

1" 

D  X> 

J  )  *  X  ( .) 

.'  x  rc  (  y  >  ai:u  J  ' 

xx<  j 

> 

c 

0 II 

it  it: 

6,100) 

IT  J  )  ,  T('( 

.1  )  ,  A 

!<  •;  .i ) 

,  ALPI-K  1  ) 

1 0  0 

r opm a r 

(  /  ,  SX  , 

t;  1. I*  f  4  x  t 

L  1  2 

,  4  .X 

, E 1 2  r,4X, 

Cl 2  C,4.X,f 

;  i.  2  ) 

c 

c 

CA 

iLL  O 

N  TOE 

PONGE  KUTTA  S 

in  or 

ION  METHOD 

TO 

c 

•j  r; 

ILVE 

IliE  GOVERNING  E 

QUA  r 

r  on 

c 

WR 

t  rr: 

(  L  ,  1  0  1 

)  X  (  1 ) 

1 1)  1. 

I  c 

IRMA  r 

in;1  s 

; ) 

II 

»:  x  c  j 

)  l  u 

ro  )  uj 

to  :: 

0 

1,3 

CAL!..  R 

UNGT  ( 

N  ,  T  ,  E  ,  X  ,  H 

!  >  M  .  K 

,  j  H 

, TEMP , AC  ,  D 

:,  Ai.ro,  TC,( 

;o ) 

N  ~  t 

GO  !  0  (14, 2 0  >  , !( 

c  c  o  v r  n  i  n  g  r:  q  u  a  r  i  o  n 

14  F'-(  (IT  J  )*r<  J)  )  AfKJ))*X<  J)  TC(J)#X(J">*X.J)  <D 

GO  ! 0  13 

P 

c  write  our  Tin:  resui  iamt  time 

c:  writs:  (c,2qo)j,x< j>  ,  t<j>  ,a< j  >  ,r 

;::o  continue 

ooo  roiiMAr  (/,sx,i3,4X,ria  C,4x,ri2  s,4X ,1:1.2. o.4x, no. s> 

c 

ir(X(J)  GT  0.00000)  GO  10  21 
c  move  id  next  temperature  step 

J = J  *  1 

r  .11  (  )  GT  100)  GO  TO  21 

A ( J ) -A <  J - 1 ) 

;<(  S  1  ■■■  X (  )  ■). ) 

1  (  j  -  r  •  j  1 ) 
i  c  ■:  .I )  =  1  c  <  j  -1 » 

V '  J  )  -B(  J  1 ) 

AB<  J >  -ALK  J-i  ) 

ALPH(  J ) -ALPH< J  1  > 

GO  TO  13 

c 

C  END  or  MARCHING  PROCESS 

2.1  CONTINUE 

1: 

C  T. P L o  r  =  1  F  OR  P L 0  TTING  TEMP.  V S  .  T I M E  G RAP H 

C  IPL0r=2  FOR  PLOTTING  CONCN .  VS.  TIME  GRAPH 

IF(IJ  EQ.I)  ;.PLOT  =  i 
rr:  U  1:0  2)  :. PLOT “2 
C  WRITE(/./i)  TO,  ALP!  I  <  1 ) 

r,  r  or: ma T  (//,  1  ox , 40 <'*'>,//,  isx , '  rn== ' , E7  3 ,4X,  '  alp i :a~ ' , r«,  .1  , / / , 

*  10X, 40 ('*'), //> 

I  i(i j  eg  1 )  wr  .1: 1  e  ( & ,  1 ) 

I I  <  r  r  eg  n  wpitc(G,2)  <  m.>  ,cocl.' ,xu_)  ,l=i, m 
j  -j  ■  1 

DO  60  L--1  ,  J 
I  1  (  L  .’  1  ( L ) 

rrcriMoi  eg  d  ti<L>-x<L> 
minor  to  2 )  ci.'L)=co<l> 

60  CONTINUE 

1.  FORMAT!//, 13X,  '  T  I  ME  (  SECONDS  )',  7X  ,'  CONCN  (CO)'  ,  10X  ,  'TEMPERATURE 

2  FORMA  T  (  1 4X  ,  E  7 . 2 , 1  0 X  ,  E.  1 2  S ,  i  0 X ,  E 1 2  .  S ) 

ir<  IPLOT  EQ.i  AND  I  t  . EQ.i)  CALL  GRAPH! J  , Ti ,Ci ,11,21 ,9  0 , S  5 , 

(?  0.0,0  0,0  0,0  0, 'TIME  IN  SECONDS > ' , ' TEMPERATURE (CENT 
@  'TEMP.  VS  TIME  GRAPH,',','  ) 
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VA  , 


[ruT'.oT  lq  2  and  ii  cu  i>  (.  all  graph <  j,ti  ,ci  ,  n  ,21  ,o  .  o  ,o .  s, 

1  0  0,0  0,0  0  ,  0  0  , 

*  M  i: ML  [ N  SECONDS,  '  ,  M  RA/  [L  COMPN  ,  '  , 

*  'CUNCN  VS  T  IflE  GRAPI  I ) 

iraiGT  i>  cau.  graphix  j  ,  u  ,ci,n,2i , '  > '  > 

LOO  CONTINUE 
1000  CONTINUE 

2  0  0  0  CONTINUE 

STOP 

END 

SUDROUT [NE  OUNCE  <N,T,r,X,i:,M,K,J,A, TEMP , AD , B , ALP H , TC , CO ) 

C  Till Ci  ROUTINE  PERFORMS  A  RUNGE-KUTTA  CALCULATION 

C  IY  GILLS  METHOD 

DIMENSION  ALPIK2G0 ) , AD (260 ) ,B<260 ) , CO (260) ,TC<260) 

D  (MENS ION  TT260) , X ( 260 ) , A < 260 ) , TEMP <260 > 

D O UDLE  PRECISION  T , F , X , T E MP , A , E , O , A D , D , A L P 1 1 , T C , C 0 
COMMON  ROI,RO,AI, DI , TR , TO , GH , V , D , Th 
N-l 
M  "M  i  1. 

GO  TO  (1,4, 0,3, 7  ) ,  M 

1  0=0  CO 0  0  0 

E  =  0  .  O 
GO  TO  9 

3  £=i  0  0  (-1  .  0  0/SQRT  (2.0) 

4  T  <  J  >  —  E  <  J  )  +  0  .  0*11 

O  X< J)=X< J) IE*(F*M  Q) 

Q=2*E#li*F  +  <  i  .  0-3 . 0*E)*Q 
E=i  .00  1  .  0 0/SQRT <2 . 0  ) 

GO  TO  9 

7  X<  J)=X(J>+H*F/£>  0/3 

TEMP ( J)=X< J)+TM 

A  (  J)~A<  i  H  RQ*CH*V#  <  TEMP  <  J)  TO  ( IK#I ( J )  )  >/  (K'J  L  ¥  A I  T  P  .1.  /  2  .  U  > 
C2=R0#SII#(TEMP<  J)  TO~<TR*T<  J)  )  >/<ROI#SI!> 

C  A<J)=A(i)/(i . 0-C2) 

C  O  <  J  )  •■=  A  <  J  )  *  D I  /  U.'Z 

C  WR I TE 6 ,  i  i  )  CO  <  J  )  ,  A  <  J  )  ,  TEMP  <  J  )  ,  II  ,  F ,  O 

1  i  FORMAT!/, SX  ,'***#**  CO='  ,  I"  12 . 4,2X,  '  A*'  ,Ei2 .4,2X, '  TEMP™ ' ,  El.  0  .  4 

*  2X ,  11=' ,T  10  4,/,5X, ' F= ' ,ri0 .4,2X, 'G=' ,Fi0 . 4,/) 

ALPIK  J)=ALP!l<i>*0  .  SK(<  1  .  0  0  0000  H  .  00 0 00 0/DEXP < 40 . #CQ < J >  )  > 

C  WRITE  (6,102)  ALPH(J) 

102  FORMAT  (1*12.0) 

D  <  J  )  =  A  L.  P 1 1  <  J  )  *  T  R  /  <  R  O I  %  A I  *  D I ) 

TC(.T)=D<  J)/D 

AD<  J)--(  ALPIK  J)/<RO*CII#V)  >#<A<  1 )  <  ROSSI  1*0#  (  TM  TO  )  /  (  ROI#AI#BI )  ) 

M  =  0 
l<  =2 

GO  TO  10 
9  K  ~  1 

.1  (J  RETURN 

END 

//GO. ET14E001  DO  UNIT=SYSSQ, 

//  D I SP ~ < , P ASS ) , DSN*44SM , SPACE* <000, < 450,50)) , 

/  /  D  C  D  "•  <  R  E  C  E  M = V  B  S ,  D  L  K  S I Z  E = 8  0  0 ,  L  R  E  C  L  =  79 6 ) 

//GO.STSIN  DD  * 

?  9 

1 2 .  10.  O .  6 .  4  2 . 

1  . 

J  90  0.70  0.60  0  5j  0.20  0.10 

0  oo 

/% 

//PLOT  EXEC  SIMP LOT 
//srsiN  DD  * 

ERAZIL  ICE  STUDY 

/* 


Figure  Cl  ( cunt  d).  Variation  of  temperature  rise  attributable  to  frazil  ice  formation  with  initial 
value  of  heat  transfer  coefficient,  aQ.  Initial  temperature  of  supercooling  is  Ta . 
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Figure  C2  (cant d).  Variation  of  temperature  rise  attributable  to  frazil  ice  formation  with  initial 
value  of  heat-transfer  coefficient,  ar).  Grid  heating  f  =  0.  Initial  temperature  of  supercooling  is  T0. 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Ettema,  R. 

Frazil  ice  formation  /  by  R.  Ettema,  M.F.  Karim  and 
J.F.  Kennedy.  Hanover,  N.H.:  U.S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory;  Springfield,  Va. : 
available  from  National  Technical  Information  Service, 
1984. 

vi,  44  p.,  illus.;  28  cm.  (  CRREL  Report  84-18.  ) 

Bibliography:  p.  28. 
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